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I.   BACKGROUND 

A.   GENERAL 

The  modern,  complex  flight  simulator  has  assumed  an 
expanded  role  in  aviation  research  and  development  as  well 
as  in  military  and  civilian  training  programs.   In  the  face 
of  the  current  energy  crisis,  manpower  and  budgetary  cut- 
backs, and  reduced  flight  operations,  the  Department  of 
Defense  (D.O.D.)  has  recognized  the  benefits  of  judiciously 
substituting  simulator  "flight"  time  for  aircraft  hours. 
Supplementing  existing  training  programs  with  instruction 
in  flight  crew  trainers  has  long  been  an  effective  means  of 
training  safe,  professional  aviators.   D.O.D.  has  establish- 
ed an  objective  of  reducing  flight  operations  25%  by  the  end 
of  Fiscal  Year  (FY)  1981.   The  Chief  of  Naval  Operations 

(CNO)  has  recently  authorized  the  substitution  of  simulator 
time  for  ten  percent  of  total  annual  flight  time  requirements 

(25  percent  of  the  annual  instrument  flight  time) .   Reference 

1  states  the  CNO's  policy: 

"As  additional  simulators  become  available  and  more 
is  learned  on  the  'transfer  of  learning'  gained  through 
the  use  of  simulators,  this  program  will  be  expanded." 

Increased  emphasis  on  the  use  of  synthetic  training 
devices  has  accelerated  the  procurement  and  funding  of 
special  part-task  trainers  such  as  the  Navy's  Air  Combat 
Maneuvering  device  and  Universal  Night  Carrier  Landing 
Simulator.   Total  D.O.D.  FY77  procurement  funding  for  train- 
ing devices  is  298.7  million  dollars  as  compared  to  88.5 
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million  in  FY74.   Research  and  development  funding  during  the 
same  period  has  nearly  doubled  [Ref.  2].   The  development  and 
deployment  of  modern  combat  aircraft  today  is  paralleled  by 
the  production  of  high-fidelity  weapon  system  trainers.   Ad- 
vanced weapon  platforms  such  as  the  F-14,  S-3A,  P3-C,  SH-2F, 
and  the  SH-3H  have  (or  will)  become  operational  concurrently 
with  advanced,  full  task,  motion  base  trainers. 

The  following  section  reviews  the  current  "state  of  the 
art"  of  the  flight  simulator  as  a  tool  in  the  training  en- 
vironment, research  and  development  community,  and  in  the 
field  of  accident  investigation.   It  references  selected 
studies  and  events  that  demonstrate  the  devices'  capabilities 
and  effectiveness. 

1.   Training 

Understanding  the  benefits  of  flight  crew  trainers  is 
essential  to  their  cost  effective  use  within  comprehensive 
training  programs.   Efficient  use  of  simulators  reduces  the 
expense  of  fuel,  weapons,  support  equipment,  and  manhours. 
It  also  reduces  the  competition  for  hardware  between  training, 
operational,  and  maintenance  departments.   Alloted  training 
hours  can  be  devoted  entirely  to  the  task  at  hand,  eliminating 
time  for  such  activities  as  pref light,  refueling,  and  clear- 
ance acceptance.   As  well  as  reducing  accidents,  it  provides 
a  safer  environment  for  practice  of  hazardous  or  emergency 
flight  operations.   The  advanced  system  incorporates  fea- 
tures such  as  "freezing"  the  problem  for  real-time  critique, 
replaying  the  maneuver  allowing  the  student  to  objectively 
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review  it  while  fresh  in  his  mind,  and  computation  of  purely 
objective  performance  measures  using  computer  based 
algorithms.   (See  Ref.  3  and  4  for  examples.) 

Instruction  in  training  devices  provides  better  control 
of  psychological  factors  such  as  stress  and  workload  levels. 
Use  of  a  select,  limited  number  of  experienced  instructors 
enhances  standardization  of  undergraduate  and  "fleet"  pilots. 

The  Naval  Air  Training  Command,  in  1972,  began  evaluation 
of  a  computer  generated  imagery  system  which  complemented  the 
existing  TA-4  2F-90  flight  simulator.   The  device  provided 
visual  scenarios  such  as  aircraft  carrier  landings,  a  bomb- 
ing range  complete  with  scoring  system,  and  carrier  catapult 
launches.   The  resultant  success  in  the  Familiarization  Stage 
of  the  advanced  student  flight  training  syllabus  was  expressed 
in  the  following  statements  [Ref.  5]: 

"The  major  area  of  benefit  has  been  the  practice 
provided  in  VFR  procedures  and  techniques.   The  addition 
of  a  carrier  visual  presentation  has  made  student  in- 
struction in  CQ    [carrier  qualification]  possible  with 
much  the  same  results  as  those  experienced  in  FAM  stage. 
The  real  benefit  of  the  simulator  however  has  proven  to 
be  in  the  weapons  delivery  phase  of  the  syllabus.   Bomb 
instructors  have  noticed  a  significant  improvement  in 
the  students'  ability  to  fly  the  pattern  and  develop  the 
principles  of  bombing  as  a  result  of  incorporating  simu- 
lator flights  into  the  weapons  stage.   ...After  a  short 
period  of  instruction,  it  was  evident  that  students  who 
had  had  several  hops  in  the  simulator  were,  to  a  signi- 
ficant degree,  outperforming  those  who  had  not." 

"The  real  key  to  the  value  of  the  simulation  is  the 
hit-spotting  program  which  provides  the  instructor  with 
real  time  readouts  of  release  altitude,  release  mach 
number,  dive  angle,  and  hit  position.   ...The  instructor 
may  also  introduce  any  of  a  multitude  of  factors  affect- 
ing mission  completion  including  adverse  weather,  air- 
craft malfunctions,  variable  winds  and  many  others 
affecting  flight  safety.   ...The  system  has  enabled  the 
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instructors  to  provide  the  student  aviator  with  a  high 
degree  of  realism  and  a  method  of  evaluation  which, 
until  now,  has  been  unavailable  even  in  the  aircraft." 

2 .   Research  and  Development 

The  use  of  cockpit  simulators  in  the  research,  devel- 
opment, test,  and  evaluation  phases  of  new  aircraft  allows 
the  aeronautical  engineer  to  develop  flight  systems  and  pro- 
cedures without  use  of  the  prototype  aircraft.   Concurrent 
development  of  mockups  provides  for  efficient  aircrew  station 
design  and  analysis  from  a  system  and  human  factors  stand- 
point.  Design  of  displays  and  controls  is  simplified  with 
construction  of  full-scale  prototype  crew  stations.   Realis- 
tic, motion  base  simulators,  driven  by  mathematical  models 
derived  from  wind  tunnel  testing,  allow  test  pilots  to  opti- 
mize aircraft  performance  capabilities  and  "handling  qualities" 
specific  to  the  mission. 

Development  of  the  U.S.  Air  Force/Boeing  Advanced  Medium 
STOL  (short  takeoff  and  landing)  Transport  was  aided 
extensively  by  the  use  of  Boeing  fixed-base  and  NASA  motion 
base  flight  simulators.   The  following  advanced  aircraft 
design  concepts  were  evaluated  on  the  devices  before 
incorporation  into  the  prototypes  [Ref .  6] . 

■  engine  bleed  air  for  leading  edge  flap  boundary 
layer  control  (BLC) 

■  redundant  mechanical  and  electrical  flight  controls 

■  aerial  delivery  procedures 

■  satisfactory  handling  qualities  and  criteria  for 
engine-out  STOL  approaches 
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■   "conventional"  piloting  techniques  for  STOL  speed 
and  flight  path  control. 
An  added  benefit,  recently  used  in  the  development  of  the 
F-14,  is  the  capability  to  design,  evaluate,  and  optimize 
specialized  subsystems.   The  Automatic  Carrier  Landing  Sys- 
tem (ACLS)  of  the  F-14  "TOMCAT"  was  developed  using  a 
3-degree  of  freedom  (pitch,  roll,  and  heave)  NASA  cockpit 
simulator  [Ref .  7] .   After  "proper  validation"  of  the  flight 
dynamics,  optimization  of  the  parameters  of  the  ACLS ' s  three 
main  components  began.   Use  of  simulated  flights  in  lieu  of 
aircraft  flights  yielded  "significant  increases  and  flexibil- 
ity in  the  number  of  combinations  [of  parameters]  that  can 
be  examined  for  aircraft  subsystems  which  have  many  parameters 
to  specify. "   Other  advantages  of  the  program  were  soon 
realized: 

"1.   Allowed  seven  iterations  of  SPN-42  [shipboard 
guidance  equipment]  parameter  values  to  be  examined  in 
two  to  three  days  instead  of  the  normal  two  to  three 
weeks  and  several  flights. 

2.  Provided  qualitative  predictions  of  pilot  com- 
ments on  closed  loop  ACLS  control  characteristics  and 
quantitative  prediction  of  performance  statistics,  pilot 
AQR  (ACLS  Quality  Rating) ,  and  to  a  lesser  extent,  open 
and  closed  loop  frequency  responses. 

3.  Saved  $120,000  and  $150,000  and  approximately 
six  months  in  the  SPN-42  optimization  phase  of  ACLS 
development . " 

"The  impact  of  this  new  T&E  capability  is  indicated 
by  the  inclusion  of  specific  periods  of  simulator  util- 
ization in  the  ACLS  development  plan  for  the  S-3A  airplane. 
Similar  utilization  will  likely  be  specified  for  the  F-18 
airplane.   In  addition,  the  same  NASA-Ames  Research  Center 
R&D  simulator  used  for  the  F-14A  has  been  reprogrammed  to 
attack  long  standing  ACLS  problems  of  the  A-7E  and  A-6 
aircraft. " 
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3.   Accident  Investigation 

Cockpit  simulators  have  proven  to  be  an  indispensible 
tool  for  reconstruction  of  pre-accident  histories.   The  pres- 
ent capability  to  use  data  from  inflight  recorders  to  derive 
flight  profiles,  environmental  conditions,  and  signals  with 
which  to  "drive"  simulators  is  impressive.   Current  digital 
recorders  can  record  up  to  64  flight  parameters  per  second. 
However  considerable  re-formatting  of  the  data  is  then  re- 
quired in  order  that  simulators  can  use  the  data  directly. 
A  current  NASA  program  is  in  pursuit  of  the  capability  to 
recreate  pre-accident  flight  conditions  using  information 
from  flight  recorders,  voice  recorder  data,  radar  tracks,  and 
meteorological  data. 

The  investigation  into  the  crash  of  an  Eastern  Air- 
lines B-72  7  in  June,  19  74,  by  the  National  Transportation 
Safety  Board  included  evaluation  of  the  wind  shear  and  down- 
draft  activity  encountered  by  the  aircraft  [Ref.  8].   Four- 
teen experienced  pilots  flew  a  total  of  54  approaches  in  a 
fixed-base  B-727  cockpit  simulator  programmed  with  four 
different  hypothesized  wind  models.   Even  with  an  a  priori 
knowledge  of  the  type  and  severity  of  the  storm  activity, 
18  flights  "crashed",  while  only  five  flights  resulted  in 
placement  of  the  aircraft  near  the  runway  threshold.   Subse- 
quently constructed  flight  profiles  (airspeed  and  altitude 
traces)  closely  resembled  those  from  the  B-727' s  flight 
recorder,  validating,  to  some  extent,  the  wind  shear  models. 
At  the  time  of  the  crash,  the  flight  crew  was  thought  to  be 
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attempting  visual  contact  with  the  runway  environment.   Seven 
of  ten  pilots  who  commented  following  the  tests  felt  that 
their  efforts  to  "go  visual"  at  decision  height  delayed  their 
recognition  of  the  sink  rate  produced  by  the  wind  shear  and 
downdraft  conditions. 

A  committee,  chaired  by  the  Director  of  Aviation 
Safety  Programs,  Naval  Postgraduate  School,  was  formed  in 
1975  to  investigate  present  use  and  potential  of  flight  sim- 
ulators in  aircraft  accident  prevention  and  investigation. 
Representatives  from  government,  NASA,  and  the  airline  indus- 
try recommended  further  research  into  the  use  of  existing 
inertial  navigation  equipment  to  record  available  flight 
accelerations  in  a  form  adaptable  to  simulators  [Ref.  9]. 
The  committee  offered  these  recommendations: 

"Crash-protected  digital  flight  recorders  should 
be  installed  in  high  valued  aircraft  (with  inertial 
navigators) .   There  should  be  a  unified  effort  to 
develop  and  implement  standardized  data  collection  and 
processing  procedures  which  would  facilitate  recreation 
of  accident/incident  conditions  in  flight  simulators. 
Voice  recorders,  especially  on  multi-crew  aircraft, 
should  also  be  considered." 

3.   NATURE  OF  THE  PROBLEM 

1.   Pilot-Simulator  Feedback  Loop 

While  accomplishing  a  certain  task  in  a  flight  sim- 
ulator, the  pilot  responds  to  various  visual,  kinesthetic, 
audio,  and  vestibular  feedback  cues  (see  Fig.  1) .   Feedback 
may  be  classified  as  either  "intrinsic"  (e.g.,  control  "feel" 
or  amplitude  of  movement  necessary  in  the  operation  of  a 
control)  or  "extrinsic"  (the  consequences  or  results  of 
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control  operation) .   The  problem  encountered  by  a  designer  is 
one  of  cost-effectiveness,  i.e.,  which  cues  are  necessary  and 
what  degree  of  fidelity  is  required  of  each  cue  to  ensure  a 
positive  transfer  of  training  at  the  particular  stage  of 
learning  for  the  given  task. 

The  requirement  for  motion  fidelity  in  a  flight  sim- 
ulator has  been  an  ongoing  question  for  years.   It  cannot  be 
debated  that  motion  cues  simulating  airframe  buffeting  as 
well  as  in  response  to  flight  control  inputs  adds  greatly  to 
psychological  fidelity,  especially  during  early  stages  of 
training.   Many  research  studies  have  compared  fixed-base 
and  motion  base  simulators  measuring  pilot  performance, 
learning  rates,  and  subjective  pilot  responses.   Where 
motion  cues  play  an  important  role  in  the  pilot's  feedback 
loop,  the  addition  of  motion  cues  creates  a  significant 
improvement  in  performance.   For  example,  during  instrument 
approaches,  the  inherent  lag  in  cockpit  instrumentation 
forces  the  pilot  to  rely  on  the  vestibular  sense  to  a  higher 
degree.   (A  possible  argument  for  fixed-base  instrument 
training,  where  the  pilot  is  forced  to  depend  on  artificial 
cockpit  references).   During  visual  carrier  landings,  the 
first  sensation  of  rapid  glide  path  deviation  is  typically 
through  the  vestibular  and  kinesthetic  sensors,  quickly 
confirmed  by  movement  of  the  visual  landing  aid.   Another 
scenario  in  which  motion  cues  are  important  to  the  pilot  is 
low-level,  high  speed  terrain  following  flights.   Where  motion 
tends  to  degrade  pilot  performance,  the  addition  of  these  cues 
allows  a  quicker  transfer  of  learning  to  the  aircraft. 
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2.   Physiology  of  Motion  Detection 

Man  perceives  the  sensation  of  motion  through  three 
main  sensors.   The  eyes  view  the  changing  world  around  him. 
The  kinesthetic  sensors  within  his  body  respond  to  the  rate 
of  change  of  movement  of  his  body  and  its  members.   The  semi- 
circular canals  and  otolith  organs  are  the  main  components 
of  the  highly  specialized  vestibular  system  of  the  inner 
ear.   The  inner  ear  canals  react  to  changes  in  rotational 
velocity,  i.e.,  angular  acceleration.   The  utricles  and 
saccules  of  the  otolith  system  respond  to  linear  accelerations 

The  present  intent  is  not  to  analyze  in  depth  the 
physiology  of  the  vestibular  system,  but  to  review  the  mo- 
tion to  which  the  organs  respond.   (A  more  comprehensive 
review  of  vestibular  physiology  and  its  function  is  available 
in  Ref.  10.) 

The  fluid-filled  semicircular  canals  form  a  roughly 
orthogonal  set  of  sensors  in  each  inner  ear.   The  endolymph 
fluid  lags  behind  the  movement  of  the  canal  walls  when  the 
head  undergoes  rotation.   This  relative  motion  tends  to  dis- 
place the  cupula  which  obstructs  an  expanded  portion  of  the 
canal  called  the  ampulla.   Deflection  of  the  cupula  excites 
sensory  hair  cells  at  its  base  causing  an  increase  in  the 
firing  rate  of  the  afferent  nerve  fibers.   Orientation  of 
the  canals  is  such  that  those  on  the  left  side  of  the  head 
are  coplanar  with  those  on  the  right  side  of  the  head, 
sensing  motion  about  the  same  axes,  but  in  the  opposite 
direction,  thereby  acting  as  a  single  unit.   It  can  be 
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assumed  that  the  central  nervous  system  responds  to  the 
difference  of  their  respective  afferent  responses. 

The  otolith  system  of  the  inner  ear  is  comprised  of 
two  utricles  and  two  saccules.   Each  is  a  mass  of  calcium 
carbonate  crystals  suspended  by  sensory  hair  cells.   Linear 
accelerations  tend  to  displace  the  mass  exciting  the  af- 
ferent nerve  endings  at  the  base  of  the  hair  cells.   The 
utricle  organs  are  oriented  so  as  to  be  sensitive  to  the 
horizontal  component  of  specific  force.   The  saccule  organs, 
approximately  perpendicular  to  the  utricles,  sense  vertical 
shear  forces  and  are  roughly  one-half  as  sensitive. 
3 .   Vestibular  Modelling  Research 

The  need  for  development  of  a  model  of  the  human 
vestibular  system  can  be  illustrated  best  when  one  considers 
the  advances  made  in  modern  aerospace  technology.   The 
high  G  environment  of  combat  aircraft  and  space  launch 
vehicles,  the  prolonged  weightlessness  of  spacecrafts,  and 
the  motion  experienced  in  large  ships  and  tall  buildings  are 
products  of  twentieth  century  technology  requiring  the 
understanding  of  man's  capacity  to  function  in  these 
environments . 

The  "state  of  the  art"  of  current  modelling  capabil- 
ities is  best  illustrated  in  a  PhD  thesis  by  Ormsby  [Ref.  11] 
At  the  MIT  Man-Vehicle  Laboratory,  under  a  NASA  grant,  Dr. 
Ormsby  developed  a  mathematical  model  based  on  the  known 
physiology  of  the  human  vestibular  apparatus.   The  signal  in 
noise  model  is  capable  of  predicting  near  threshold  and 
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supra-threshold  perceived  motion  and  attitude  in  the  pres- 
ence of  random  noise  and  spontaneous  firing  of  the  afferent 
sensors.   It  was  then  refined  to  be  consistent  with  current 
neurophysical  and  psychophysical  data.   Included  is  a  model 
of  the  central  processor  (the  brain)  and  its  interpretation 
of  afferent  dynamics  with  an  a  priori  knowledge  of  the  mo- 
tion stimulus.   Since  it  was  difficult  to  obtain  data 
resulting  only  from  vestibular  cues,  the  model  is  said  to 
also  reflect  responses  resulting  from  tactile  and  kinesthetic 
motion  cues. 

Separate  models  of  the  semicircular  canals  and  the 
otolith  organs  are  combined  to  yield  the  following  estimates 
in  the  coordinate  system  of  the  head: 

■  perceived  linear  acceleration 

■  perceived  rotational  velocity 

■  perceived  attitude  (a  unit  vector  in  the  perceived 
"down"  direction) . 

The  model  has  been  successful  in  predicting  illusions 
such  as  the  dynamic  elevator  illusion  in  the  absence  of  vi- 
sual inputs.   It  also  has  produced  results  in  close  agreement 
with  data  describing  perceived  lateral  tilt  angle  as  a  func- 
tion of  actual  tilt  angle  in  "1G"  and  "2G"  environments.   The 
pitch  and  roll  sensations  during  a  catapult  launch  were 
investigated  and  used  to  further  refine  the  model. 

It  must  be  emphasized  that  the  version  of  the  model, 
as  used  in  this  thesis,  is  strictly  a  vestibular  model  with- 
out visual,  tactile,  or  kinesthetic  inputs.  The  integration 
of  vestibular  and  visual  cues  is  the  subject  of  an  ongoing 
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research  program  at  the  MIT  Man-Vehicle  Laboratory 

[Refs.  12  and  13] . 

4 .   Components  of  a  Motion  Base  Control  System 

The  human  vestibular  system  is  sensitive,  not  to 
motion,  but  to  the  higher  derivatives  of  motion  such  as 
linear  and  angular  acceleration.   Thus  only  the  "rate  of 
change"  of  motion  need  be  simulated.   The  simulation  prob- 
lem, then,  is  how  to  optimally  simulate  acceleration  cues 
given  fixed  operating  limits.   (Hardware  has  inherent  size, 
velocity,  and  inertial  limitations,  while  computers  have 
speed  and  accuracy  limits.)   A  computer  model  of  the  air- 
craft produces  acceleration  vectors  in  response  to  control 
inputs  and  external  factors  such  as  wind  and  turbulence 
which  are  converted  to  simulator  drive  commands.   Linear 
accelerations  are  simulated  by  tilting  the  cockpit  cab 

(such  as  backward  on  application  of  thrust)  at  a  rate  be- 
low human  rotational  thresholds.   Thus  lateral  and  longitu- 
dinal acceleration  cues  can  be  simulated,  but  only  the 
higher  frequency  vertical  cues  (such  as  turbulence  and 
rapid  pullups)  can  be  simulated.   After  the  cue  has  been 
transmitted,  the  cab  is  gradually  returned  to  its  neutral 
position.   This  tendency,  known  as  "wash-out",  keeps  the 
cab  near  its  center  position  to  maximize  its  operating 
envelope  for  subsequent  cues. 

This  section  describes  the  main  components  of  motion 

base  control  systems.   The  vast  majority  of  motion  base 

simulators  have  either  "synergistic"  platforms  supported  and 
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driven  by  linear  hydraulic  actuators  or  electro-mechanical 
drive  systems  with  independent  degrees  of  freedom.   Digital 
computer  programs  for  both  system  types  are  comprised  of 
three  main  components : 

■  aircraft  models  which  produce  the  particular 
flight  dynamics  and  control  feedback  ("feel") 
characteristics 

■  "wash-out"  filters  which  remove  low  frequency 
accelerations  which  tend  to  generate  large 
amplitude  movements 

■  "drive"  programs  which  transform  the  aircraft 
motion  to  simulator  motion  commands . 

This  sequential  process  is  shown  in  Figure  1.   The  motion 
of  the  aircraft  about  its  six  degrees  of  freedom  is  used  to 
drive  cockpit  instrumentation,  update  navigational  equipment, 
and  drive  external  visual  displays.   The  drive  program 
consists  of  algorithms  which  compute: 

■  residual  tilt  which  uses  gravitational  forces 
to  simulate  linear  accelerations 

■  lead  compensation  for  "quickening"  the  position 
commands  to  compensate  for  inertial  lag  in  the 
machine's  dynamics 

■  cab-to-drive  mechanism  transformations 

■  "soft"  position,  velocity,  and  acceleration 
limiting  networks. 

Although  many  additional  safety  features  and  utility  func- 
tions are  present  in  most  programs,  these  are  the  primary 
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motion  generating  components.   (See  Refs.  14  and  15  for  more 
in-depth,  theoretical  descriptions  of  control  systems.) 
5.   Simulator  Fidelity 

Studies  which  evaluate  motion  control  systems  employ 
sophisticated  means  of  calculating  performance  measures  and 
compiling  subjective  pilot  responses.   Numerous  studies  have 
evaluated  the  handling  quality  differences  between  fixed 
and  motion  base  simulators  in  an  effort  to  understand  the 
need  for  and  the  contribution  of  motion  cues  to  the  pilot. 
(See  Refs.  15  and  16  for  examples.)   Improving  stability 
characteristics  for  new  aircraft  or  simulators  often  in- 
volves the  analysis  of  pilot  performance  while  system  param- 
eters are  systematically  changed.   Verbal  responses  are 
analyzed  to  discern  which  aspects  of  the  motion  cues  are 
deficient.   Design  of  an  aircraft  control  system  typically 
entails  a  series  of  compromises  and  tradeoffs  yielding 
handling  qualities  seldom  optimum  for  any  specific  operating 
regime.   Therefore,  optimum  performance  does  not  reflect  the 
fidelity  of  a  motion  system.   Approximate  duplication  of 
such  fundemental  maneuvers  as  coordinated  turn  entries  has 
been  accomplished  by  matching  aircraft  and  simulator  accel- 
erations about  all  six  degrees  of  freedom.   Since  a  precise 
duplication  is  unattainable  due  to  the  simulator's  physical 
constraints,  only  the  human  operator  can  properly  evaluate 
the  differences. 

Optimization  of  the  control  system  for  a  new  simu- 
lator is  typically  accomplished  through  numerous  flights  by 
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an  experienced  test  pilot  while  engineers  vary  stability 
parameters.   Many  unpublished  fidelity  programs  involve 
hundreds  of  manhours  and  expensive  computer  time  in  search 
of  optimal  control  parameters.   Often  pilots  need  refresher 
flights  in  the  aircraft  to  regain  familiarity  with  its 
stability  characteristics.   One  such  study  [Ref.  22]  esti- 
mated that  "eight  hours  of  simulator  time  is  sufficient  to 
destroy  a  pilot's  familiarity  with  the  airplane  to  the  point 
that  he  can  no  longer  make  valid  qualitative  judgements  on 
flying  qualities  during  closed  loop  tests." 

C.   PRESENT  PURPOSE 

Present  D.O.D.  emphasis  is  on  new,  higher  fidelity 
flight  simulators.   A  second  avenue  of  approach  is  the 
improvement  of  present  training  devices  having  less  than 
optimal  flying  qualities.   Lack  of  fidelity  has  reduced  the 
purpose  of  many  devices  to  mere  procedural  trainers,  while 
the  motion  and  visual  systems  go  unused.   With  new  technology 
visual  display  systems,  the  shortcomings  of  the  handling 
qualities  and  performance  characteristics  are  magnified. 
As  simulators  assume  a  greater  role  in  flight  training  and 
proficiency  programs  teaching  advanced  flying  skills,  the 
need  for  realistic  duplication  of  the  aircraft  is  apparent. 

The  Naval  Air  Test  Center  (NATC)  has  been  conducting  a 
fidelity  improvement  program  for  the  Navy's  inventory  of 
Operational  Flight  Trainers  and  Weapon  System  Trainers 
[Ref.  22] •   A  team  of  fleet  experienced  pilots,  flight  test 
engineers,  simulator  and  computer  experts  have  been 
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successful  in  identifying  performance,  stability,  and  con- 
trol parameters  contributing  to  poor  fidelity  in  all  the 
subsystems.   Although  some  hardware  changes  have  been  made, 
the  majority  of  fidelity  improvements  have  been  accomplished 
by  software  reprogramming. 

"Aircraft  manufacturers  generally  rely  on  wind 
tunnel  data  as  a  data  base.   Flight  test  data  are  not 
available  because  the  first  training  device  is  usually 
delivered  to  the  Fleet  at  about  the  same  time  that  the 
first  airplane  comes  off  the  production  line.   Since 
wind  tunnel  data  represent  at  best  only  an  estimate  of 
airplane  flight  characteristics,  the  use  of  these  data 
as  a  data  base  for  a  simulator  results  in  poor  flying 
qualities  simulation.   In  addition,  there  has  been  no 
organized  effort  within  the  Navy  or  by  any  manufacturer 
to  reprogram  a  simulator  once  flight  test  data  are 
available . " 

"Recent  advances  in  the  practical  applications  of 
control  and  estimation  theory  have  made  it  possible  to 
extract  stability  derivatives  and  other  parameters 
from  flight  test  data  to  a  degree  and  with  accuracy 
not  possible  in  the  past." 

The  degree  of  realism  required  in  equipment,  visual  dis- 
plays, and  motion  systems  is  a  function  of  the  level  of 
training  being  performed  and  the  amount  of  learning  transfer 
desired.   For  a  specific  training  device,  whatever  cost- 
effective  measures  (up  to  the  point  of  diminishing  returns) 
increase  its  fidelity  would  be  welcomed,  thus  improving  the 
transfer  of  learning. 

The  goal  of  this  thesis  is  the  development  of  a  technique 
which  would  allow  the  simulator  programmer  (such  as  those  of 
NATC's  fidelity  improvement  program)  to  optimize  the  archi- 
tecture and  choice  of  parameters  for  simulator  software. 
Work  on  an   ongoing  research  program  by  Dr.  L.  R.  Young  at 
the  NASA  Langley  Research  Center  and  the  MIT   Man-Vehicle 
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Laboratory  on  the  integration  of  visual  and  motion  cues 
provided  many  original  ideas  for  this  work.   An  experiment 
was  performed  at  the  NASA  Ames  Research  Center  in  the  Flight 
Simulation  Laboratory  on  a  six-degree-of-freedom  research 
simulator.   The  experiment  was  designed  to  test  the  hypothe- 
sis that  there  was  no  difference  between  simulator  motion 
cues  presented  to  the  pilots.   Five  pilots  performed  a  total 
of  161  simulator  "runs",  providing  subjective  ratings  of 
different  degrees  of  fidelity. 

Using  Ormsby ' s  model  of  human  dynamic  orientation,  a 
computer  program  was  written  to  simulate  the  pilot's  subjec- 
tive evaluation  of  fidelity  in  objective  terms.   Validation 
of  the  simulation  was  performed  using  the  experimental  re- 
sults as  a  data  base.   Another  aim  of  this  report  is  to 
confirm  Dr.  Ormsby ' s  model  to  whatever  extent  possible. 
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II.   METHOD 

A.   EXPERIMENT 
1 .   Design 

The  pilot-subjects  were  asked  to  provide  responses 
to  different  levels  of  motion  fidelity.   Fifteen  sets  of 
"wash-out"  filter  parameters  were  selected  in  an  effort  to 
produce  discernably  different  motion  cues  during  a  flight 
task.   For  each  "flight",  subjects  were  asked  to  rate  the 
rotational  and  linear  acceleration  cues  relative  to  a 
rating  scale  and  to  provide  a  combined  fidelity  score. 
Thus,  the  dependent  variables  for  the  experiment  were  the 
three  subjective  ratings  of  the  motion  cues.   The  indepen- 
dent variables  in  the  analysis  were  the  motion  cues  and 
parameter  sets. 

The  experimental  results  were  analyzed  according  to 
a  three-way  factorial  analysis  of  variance.   Two  fixed  fac- 
tors were  motion  (three  levels)  and  parameter  sets  (15 
levels).   The  third  was  subjects,  a  random  factor  (five 
levels) .   A  conceptual  model  of  the  experimental  design  is 
shown  in  Figure  2. 
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Figure  2:   Model  of  Experimental  Design 

2 .   Subjects 

Since  the  simulator  was  programmed  as  an  F-15  (a 
combat  fighter  aircraft) ,  subjects  for  the  study  were  chosen 
for  their  experience  as  pilots  in  high  performance  aircraft. 
Although  the  subjects  obtained  for  the  experiment  were  ex- 
perienced in  military  combat  aircraft,  only  one  was  relative- 
ly current.   The  other  four  subjects  were  military  trained 
airline  pilots  or  flight  engineers.   Their  pilot  time  exper- 
ience is  shown  in  Table  I.   Each  was  a  volunteer  and  was 
compensated  for   his   participation. 
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SUBJECT  HELP  SINGLE -ENGINE  MULTI-ENGINE  HIGH  PERFORMANCE  JET 

1  -  550  2400  3050 

2  -  180  860  130 

3  600  160  350  250 
4-30  -  800 
5      -  1600  7500  2150 

Table  I:  Subject  Experience  (pilot  hours) 

3 .   Stimuli  and  Apparatus 
a.   Simulator 

The  experiment  was  performed  at  the  NASA  Ames 
Research  Center,  in  the  Flight  Simulation  Laboratory.   The 
simulator  used  in  the  experiment  was  the  S.01  All-Axis 
Motion  Generator,  manufactured  by  the  Northrop  Corporation. 
The  six-degree-of-f reedom  research  simulator  has  independent 
drive  mechanisms  with  a  single-seat  cockpit  in  the  center  of 
its  operating  envelope.   The  machine's  linear  position  limits 
are  ±9  feet;  the  rotational  limits  are  ±4  5  degrees  from  its 
neutral  position.   The  simulator  was  controlled  by  an  EAI 
84  00  computer.   The  computer  programming  included  "soft" 
position,  velocity,  and  acceleration  limiting  networks  which 
prevent  the  simulator  from  reaching  the  "hard"  mechanical 
limits.   A  more  in-depth  study  of  the  machine's  operating 
characteristics  and  mathematical  development  of  its  "drive 
logic"  is  contained  in  Reference  15. 

As  with  most  simulators,  the  cockpit  was  en- 
closed preventing  visual  reference  to  outside  of  the  cab. 
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The  cockpit  configuration  was  that  of  a  typical  single-seat 
aircraft,  controlled  by  a  stick,  rudder  pedals,  and  power 
lever.   Since  visual  cues  from  external  displays  and  arti- 
ficial attitude  instrumentation  tend  to  interact  with 
vestibular  cues  in  a  way  not  fully  understood,  use  of  these 
displays  was  minimized.   (The  simulator  was  equipped  with  a 
terrain  modelboard  visual  system,  but  the  projection  equip- 
ment was  off  for  the  duration  of  the  experiment.)   The  main 
instrument  panel  was  organized  in  a  conventional  "T"  arrange- 
ment as  shown  in  Figure  3.   Although  the  basic  instrument 
array   included  an  artificial  horizon  (attitude  gyro) ,  it 
was  felt  that  its  absence  would  impose  an  unusually  high 
workload  level  on  the  subject  and  tend  to  lengthen  the  learn- 
ing  process  inherent  in  the  task.   Also,  the  more  familiar 
the  cockpit  environment  was  to  the  subject,  the  quicker 
would  be  his  adaptation  to  the  task,  allowing  maximum 
attention  to  the  motion  cues. 

An  intercom  system,  with  a  "hot"  microphone  for 
the  pilot,  allowed  two-way  communication  between  the  subject, 
system  operators,  and  experimenter.   Verbal  responses  by  the 
pilots  were  recorded  on  a  cassette  tape  recorder.   The  ex- 
perimenter monitored  the  progress  of  the  task  on  a  strip 
chart  recorder,  displaying  altitude,  airspeed,  vertical 
speed,  and  heading. 
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1  -  Clock  with  Sweephand 

2  -  Airspeed  Indicator 

3  -  Attitude  Gyro 

4  -  Altimeter 

5  -  RPM  Gauge 


6  -  Angle  of  Attack  Indicator 

7  -  Horizontal  Situation 

Indicator 

8  -  Vertical  Speed  Indicator 

9  -  Turn  Needle  and  Ball 


Figure  3:   Main  Instrument  Panel 

The  "wash-out"  filter,  which  interfaces  the 
mathematical  aircraft  model  with  the  "drive"  program, 
is  composed  of  six  separate  digital  filters,  one  for  each 
degree  of  freedom  or  "channel".   Each  second  order,  high 
pass  filter  is  described  by  a  linear  gain,  a  natural 
frequency,  and  a  damping  ratio  with  a  Laplace  transform 
equation  of  the  form: 
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Ks2 

s2  +  2£wns  +  w2 

where:   K  =  linear  gain 

s  =  Laplace  complex  variable 

£  =  damping  ratio 
wn  =  natural  frequency 

Through  consultation  with  NASA  personnel,  15 
sets  of  parameters  (18  parameters  each)  were  chosen  to 
produce  noticeably  different  levels  of  realism.   Values  of 
the  filter  parameters  are  presented  in  Appendix  B.   The 
fidelity  associated  with  each  set  could  not  be  predeter- 
mined, except  that  the  range  over  which  each  parameter 
varied  was  thought  to  be  sufficient  to  produce  noticeable 
differences  in  motion  cues. 
b.   Task 

Initially,  the  experiment  was  designed  to  ana- 
lyze two  different  tasks  in  order  to  examine  any  significant 
difference  in  subject  responses.   Due  to  time  constraints 
and  equipment  malfunctions,  the  second  task,  an  ILS  (instru- 
ment landing  system)  approach  to  a  waveoff,  was  eliminated 
from  the  study. 

Design  of  the  task  reflected  many  considerations 
It  was  recognized  that  the  task  should  be: 

■  similar  to  those  used  in  RDT&E  and  training  programs 
so  that  the  maximum  extrapolation  of  results  could  be 
made 

■  relatively  familiar  to  the  subjects  and  well  within 
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their  capabilities  to  lessen  the  time  spent  mastering 
it 

■  well  within  the  capability  of  the  simulator  to  mini- 
mize machine  error  and  breakdowns 

■  long  enough  to  expose  the  subject  to  a  full  variety 

of  motion  cues  and  allow  him  ample  time  to  subjectively 

evaluate  them. 

A  three  and  one-half  minute  "basic  airwork"  task 
was  chosen  to  simulate  the  low  frequency  flight  maneuvers 
common  to  the  terminal  phase  of  flight  operations  (climb-out 
after  takeoff,  descent  prior  to  approach  and  landing,  etc.). 
Subjects  were  provided  enough  information  to  rapidly  learn 
the  task  and  to  execute  it  consistently  from  trial  to  trial. 
As  shown  in  Figure  4,  two  versions  of  the  flight  pattern 
were  devised.   The  direction  of  the  initial  turn  was  random- 
ized.  During  the  second  replication  of  each  parameter  set, 
the  turn  was  made  in  the  opposite  direction, 
c.   Rating  Scales 

Subjects  were  asked  to  subjectively  judge  the 
fidelity  of  the  motion  cues  encountered  during  the  task  us- 
ing two  rating  scales.   The  pilots  scored  separately  the 
linear  (or  translational)  and  the  rotational  motion  cues  on 
a  graphic  eight-point  scale.   The  cues  were  also  rated  in 
comparison  to  a  "standard"  motion  (presented  using  the 
established  set  of  wash-out  parameters  employed  on  the  simu- 
lator) .   Thus,  the  subjects  rated  the  motion  fidelity  on 
both  an  absolute  and  comparative  basis. 
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4 .   Procedure 

After  the  subjects  read  the  instructions  (see  Appen- 
dix A) ,  further  questions  were  answered.   The  intent  and 
goals  of  the  experiment  were  discussed  freely  with  the  sub- 
jects.  They  were  allowed  to  become  familiar  with  the  simu- 
lator cab  during  one  practice  run  before  the  trial  runs 
began.   Each  was  given  a  pilot  kneeboard  containing  a  dia- 
gram of  the  task,  the  two  rating  scales,  and  a  list  of  topics 
from  which  relevent  comments  were  volunteered  following  each 
run. 

The  sequence  of  wash-out  parameter  sets  presented  to 
the  subjects  was  randomized;  a  different  sequence  was  pre- 
sented to  each  subject.   The  "standard"  motion  was  presented 
once  at  the  beginning  of  each  session  (twice  at  the  start  of 
a  subject's  first  session)  and  every  fourth  run  thereafter. 
A  typical  session  of  12  simulator  runs  is  shown  in  Figure  7. 

P0   P2   P8   p15   P0   Pll   P14   P6   P0   Pl   P5   P13 
where:   PQ  =  "standard"  parameter  set 

Figure  7:   Typical  Sequence  of  Parameter  Sets 

No  reference  was  made  to  the  parameter  set  number  or  their 
values  during  the  trial  runs. 

At  the  end  of  each  task,  subjects  rated  the  motion 
cues  and  provided  subjective  comments  relevent  to  the  spe- 
cific run  via  the  intercom.   The  parameter  set  was  changed 
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and  the  simulator  was  reinitialized  for  the  next  run.   After 
each  session,  subjects  were  allowed  to  review  their  responses 
and  discuss  the  session  (comments  were  not  recorded) . 

Factors  which  were  uncontrollable  were  the  number  of 
sessions  and  runs  per  subject.   Due  to  equipment  malfunctions 
and  scheduling  constraints  the  sessions  and  runs  varied 
widely  as  shown  in  Table  II. 


SUBJECT 

NO. 

SESSIONS 

NO 

.  RUNS* 

1 

3 

24 

2 

3 

37 

3 

4 

42 

4 

2 

26 

5 

2 

31 

♦Includes  "Standard  Motion"  Runs 

Table  II:   Sessions  and  Runs  by  Subjects 

B.   SIMULATION 

1.   Model  Assumptions 

The  subjective  process  by  which  a  pilot  compares  the 
motion  of  a  simulator  with  that  of  an  aircraft  is,  obviously, 
a  complex  psychological  one.   If  a  model  of  this  human  com- 
parative method  were  constructed,  it  could  predict  a  pilot's 
response  to  the  fidelity  of  a  simulator's  motion  character- 
istics.  The  structure  of  the  computer  simulation  described 
in  this  section  is  based  on  the  assumption  (as  proposed  by 
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Young,  Ref.  13)  that  the  pilot's  evaluation  is  a  function  of 
the  difference  between: 

■  the  simulator  motion  cues  as  perceived  by  the  pilot 
and 

■  the  perceived  motion  of  the  aircraft  as  experienced 
by  the  pilot. 

Assuming  that  Ormsby's  vestibular  model  of  human  dy- 
namic orientation  provides  a  reliable  estimate  of  motion 
perception,  it  was  then  a  matter  of  deriving  motion  histories 
of  the  aircraft  and  the  simulator  performing  "identical" 
maneuvers . 

2.   Design 

During  the  experiment  described  in  Section  A,  the 
output  of  the  computerized  aircraft  model  (see  Fig.  1)  was 
recorded  on  magnetic  tape  during  each   run.    The  data  con- 
sisted of  accelerations  for  each  degree  of  freedom  recorded 
during  each  cycle  of  the  discrete  programming.   The  cycle 
time  used  in  the  experiment  was  50  milliseconds,  chosen  so  as 
to  be  compatible  with  the  input  sampling  rate  of  the  dynamic 
orientation  program.   (Sufficient  data  was  also  recorded  to 
document  the  specific  pilot  and  run  number.) 

The  flight  acceleration  data  was  used  as  an  input  to 
the  S.01  simulator's  drive  program,  a  copy  of  which  is  used 
in  this  simulation.   The  simulation  is  a  deterministic  model 
with  no  random  variables.   Figure  8  is  a  block  diagram  of  the 
simulation.   Since  the  dynamic  response  of  the  motion  platform 
contributes  to  the  motion  cues  presented  to  the  pilot,  but  is 
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not  a  factor  in  this  simulation,  the  lead  compensation  net- 
work in  the  program  is  bypassed.   The  result  is  the  closest 
approximation  available  to  the  actual  motion  of  the  cab. 
The  outputs  from  this  program  are  linear  acceleration  com- 
mands and  rotational  position  commands  of  the  platform's 
drive  mechanism.   A  transformation  then  converts  these  sig- 
nals back  to  vectors  describing  the  cockpit  motion.   The 
rotational  position  commands  are  then  differentiated  to 
produce  rotational  velocity  vectors. 

Inputs  to  the  human  dynamic  orientation  program  are 
discrete  time  histories  of  specific  force  vectors  (in  units 
of  "G's")  and  angular  velocity  vectors  (radians/sec.)  in  the 
coordinate  system  of  the  pilot's  head.   It  is  sampled  at  dis- 
crete intervals  of  10  0  milliseconds;  thus  the  simulator  and 
aircraft  motion  vectors  are  sampled  at  every  other  time 
frame.   Although  this  represents  a  loss  of  input  information, 
considerable  recomputation  of  many  model  parameters  would 
have  been  necessary  to  alter  the  input  or  output  rates.   The 
output  of  Ormsby's  model  consists  of  vectors  in  the  coordi- 
nate system  of  the  subject  describing  his  perception  of 
rotational  velocity,  linear  acceleration,  and  attitude  (per- 
ceived vertical  or  "down"  direction).   The  model's  output 
rate  is  one  sample  per  second. 

Motion  vectors  of  the  aircraft  are  computed  from  the 
flight  data  after  the  rotational  velocity  is  formed  by  inte- 
grating the  rotational  acceleration.  These  vectors  are  then 
used  as  inputs  to  the  dynamic  orientation  program  to  produce 
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perceived  aircraft  motion.   Notation  used  in  the  block 
diagram  is  defined  as  follows: 


X ,  ,  X   =  aircraft  model  linear  and  rotational  acceleration 
vectors,  respectively 

SFa,  SFq  =  specific  force  vectors  for  the  aircraft  and 
simulator,  respectively 

w^,  Wg  =  rotational  velocity  vectors  for  the  aircraft 
and  simulator,  respectively 


Aa,  Ag  =  perceived  linear  acceleration  vectors  for  the 
aircraft  and  simulator,  respectively 


wHa'  ^s  =  Perceived  rotational  velocity  vectors  for  the 
aircraft  and  simulator,  respectively 

DOWN -j ,  DOWN   =  perceived  attitude  vectors  for  the  air- 
craft and  simulator,  respectively. 
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III.   RESULTS 

A.   EXPERIMENTAL  RESULTS 

1 .   Subjective  Responses 

Verbal  responses  recorded  at  the  end  of  each  task 

allowed  the  experimenter  to  understand  such  considerations 

as  the  learning  curve,  workload,  etc .,  associated  with  the 

task.   The  act  of  commenting  on  each  run  gave  subjects  the 

opportunity  to  critically  review  the  flight  and  to  adjust 

their  scores  as  opinions  were  formulated.   This  section  is 

a  brief  summary  of  pilot  comments  pertaining  to  the 

experimental  flights. 

Differentiation  between  the  rotational  and  transla- 

tional  motion  cues  was  discussed  with  the  subjects  as  the 

experiment  progressed  in  hopes  that  more  reliable  responses 

would  result.   Following  his  fifth  run,  one  subject  attempted 

to  further  define  his  task: 

"I'm  to  separate  pitch  and  roll  from  sideways,  hori- 
zontal, or  longitudinal  axis  motion  as  linear,  whether 
it  be  throttle  acceleration  or  just  sideways  gusts." 

Several  pilots  commented  that  they  were  more  confi- 
dent of  their  responses  during  second  and  subsequent  ses- 
sions as  they  familiarized  themselves  with  the  task  and 
learned  more  about  the  characteristics  of  the  motion  to  be 
rated.   One  pilot  mentioned  that  it  was  difficult  to  dis- 
cern whether  motion  cues  were  more  realistic  or  his  ability 
to  control  the  simulator  had  improved.   During  a  subject's 
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second  session  he  commented: 

....the  oscillations  I  was  in  the  other  day;  at 
first  I  thought  I  was  experiencing  pitch  oscillations. 
What  they  are  is  the  linear,  up  and  down,  that  I  was 
feeling  before;  so  I'm  getting  used  to  the  motion.... 
at  least  understanding  it  better." 

All  the  pilots  noted  that  the  workload  level  varied 
significantly  from  run  to  run  as  the  motion  cues  changed. 
Abnormal  cues  often  caused  subjects  to  react  improperly;  a 
response  which  presumably  would  not  have  occurred  in  a  fixed- 
base  simulator.   One  subject  noted  the  difference  between 
the  motion  and  cockpit  instrument  indications: 

"I  was  getting  a  definite  rotational  pitch  sensa- 
tion and  frequently  corrected  with  [forward]  stick.... 
which  I  didn't  want  to  do." 

All  subjects  noted  vibrations  induced  by  the  motion 
base  mechanism.   Normally  turbulence  masks  these  anomalies  to 
some  extent.   Since  it  was  not  used  in  the  experiment,  the 
mechanical  limitations  tended  to  put  a  ceiling  on  the  subjec- 
tive motion  ratings.   One  subject  was  allowed  a  practice  run 
with  a  minimal  level  of  turbulence  at  the  beginning  of  his 
second  session.   He  commented: 

".. .turbulence. . .simply  detracts  from  the  ability 
to  dif ferentiate. . .how  the  controls  are  functioning. 
I  think  I'm  less  able  to  determine  how  much  of  it  is 
...turbulence  or  control  response  with  the  turbulence." 

2 .   Correlation  of  Ratings 

Scores  from  the  three  rating  scales  represent  discrete 

data  possessing  ordinal  properties.   It  can  be  argued  that 

each  rating  scale  attempts  to  define  levels  separated  by  equal 

intervals;  it  is  not  immediately  evident  that  each  subject 
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perceived  the  scales  in  this  manner.   Thus,  it  cannot  be 
assumed  that  the  data  comes  from  an  interval  scale,  an 
assumption  for  use  of  parametric  statistics.   As  a  measure 
of  the  reliability  of  the  data,  nonparametric  correlations 
were  calculated.   Scores  from  the  standard  motion  were  ex- 
cluded from  this  and  subsequent  analysis  since  subjects  were 
told  which  runs  were  the  standard  and  instructed  to  use  it 
as  a  comparative  reference  for  the  other  15  parameter  sets. 
An  average  score  was  formed  for  each  subject/motion/param- 
eter set  combination  analogous  to  each  "cell"  of  the  experi- 
mental design  (see  Fig.  2) .   Spearman  rank  correlation  coef- 
ficients between  subjects  were  computed  for  the  linear, 
rotational,  and  combined  ratings  separately.   The  correla- 
tions with  respective  significance  levels  are  shown  in 
Table  III. 

SUBJECTS 


S    1     1.0 

U 

B    2 

J 

E    3 

C 

T    _4 

S 


Significance  Levels:   *p<.05,  **p<.01,  ***p<.005 


Table  III  (a);   Correlation  Between  Subjects  for 

Linear  Motion  Ratings 


2 

3 

4 

5 

5240* 

.0731 

.6810*** 

.1264 

1.0 

.1435 

.5668** 

.1722 

1.0 

.2263 

.5666** 

1.0 

.0250 
1.0 
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SUBJECTS 
12  3  4  5 

S    1     1.0       .4857*     .1244      .6081**    .4367 

U 

B    1  1.0       .2823      .5957**    .4157 

J 

E    2  1.0      -.0703      .3963 

C 

T    i  1.0       .2485 

S 

5  1.0 


Table  111(b):   Correlation  Between  Subjects  for 

Rotational  Motion  Ratings 


SUBJECTS 

12  3  4  5 

S    1     1.0       .5449*     .5065*     .7907***   .2456 

U 

B    2  1.0       .5321*     .6147**    .4700* 

J 

E    3  1.0       .4481*     .4968* 

C 

T    4  1.0       .2564 

S 

5  1.0 


Significance  Levels:   *p<.05/  **p<.01,  ***p<.005 


Table  III(c):   Correlation  Between  Subjects  for 

Combined  Motion  Ratings 


Overall  agreement  between  the  five  subjects  was  com- 
puted using  the  Kendall  coefficient  of  concordance  as  shown 
in  Table  IV. 
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RATINGS 
Linear      Rotational      Combined 
.4483*         .4835*  .5914* 

Significance  Level:   *p<.005 


Table  IV:   Coefficient  of  Concordance  Among  Subjects 
for  Motion  Ratings 


The  ability  of  subjects  to  repeat  ratings  given  the 
same  stimuli  was  investigated.   Due  to  time  constraints  only 
41  subject/parameter  set  combinations  were  replicated.   Cor- 
relations between  the  first  and  second  observations  for  these 
combinations  is  shown  in  Table  V. 

RATINGS 
Linear      Rotational      Combined 
.4882*         .4234*  .5559** 

Significance  Levels:  *p<.005,  **p<.  0005 

Table  V:   Spearman  Rank  Correlation  Between  Replications 

3 .   Analysis  of  Variance 

Analysis  of  Variance  techniques  were  used  to  deter- 
mine if  different  motion  cues  significantly  affected  the 
ratings  and  if  interactions  existed  which  would  help  explain 
the  weak  between-subject  correlations  obtained.   Although 
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the  discrete  rating  scale  data  violates  the  assumptions  of 
interval  data  and  normality  of  the  error  distribution,  an 
approximate  ANOVA  was  computed.   Since  the  rating  scale  con- 
sists of  essentially  15  levels  separated  by  intervals  de- 
signed to  be  equal,  the  assumption  of  interval  data  is  not 
entirely  untenable.   Madill  [Ref.  17]  reports  that  the  dis- 
tribution of  the  error  term  in  the  regression  model  for 
discrete  rating  scale  data  is  often  binomial  in  nature. 
Thus,  the  variance  is  not  independent  of  the  mean  as  re- 
quired for  ANOVA  analysis.   The  transformation  which  attempts 
to  equalize  the  error  variance  for  the  binomial  distribution 
is  the  arcsinyX   transformation.   A  transformation  of  the 
form: 


R'  =  7arcsinV  ~-  -  1.0  (D 

where:   R  =  raw  rating  scale  score 
R'  =  transformed  score 

was  used.   The  effect  is  minimal  in  the  range  3.0  to  6.0  and 
tends  to  expand  the  ends  of  the  scale.   The  ANOVA  table  is 
shown  in  Table  VI.   Due  to  the  large  percentage  of  missing 
data,  the  three-way  interaction  sum  of  squares  is  included 
in  the  error  sum  of  squares  with  its  degrees  of  freedom. 


49 


SOURCE 


d.f . 


3.S. 


M.S. 


F  RATIO 


Subjects  (S)  4  83.91  20.98  44.24* 

Motion  Cues  (M)  2  1.08  .54  .33 

Parameter  Sets  (P)  14  119.99  8.57  4.86* 

S  x  M  8  13.13  1.64  3.46* 

S  x  P  56  98.71  1.76  3.72* 

M  x  P  28  13.25  .47  .99 

Error  235  111.91  .48 


Total 


347    441.98 


Significance  level:   *p<.  001 

Table  VI :   ANOVA  Table  for  the  Experiment 

A  nonparametric  analysis  of  variance  was  performed 
to  support  the  results  of  the  approximate  ANOVA  above.   The 
subject/parameter  set  averages  were  tested  using  the  Friedman 
two-way  analysis  of  variance  by  ranks  [Ref.  18].   The  three 
motion  cues  were  investigated  separately.   Results  showed 
that  both  the  linear  and  rotational  motion  ratings  were  sig- 
nificantly different  at  the  .01  level.   The  combined  ratings 
were  different  at  the  .001  significance  level. 

A  Duncan  multiple  range  test  with  unequal  replications 
[Refs.  19  and  20]  was  performed  on  the  means  of  the  parameter 
sets.   Results  using  the  .05  level  of  significance  are  shown 
in  Figure  9.   Any  two  means  underscored  by  the  same  line  are 
not  significantly  different. 
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Figure  9 :   Parameter  Set  Averages 

B.   SIMULATION  RESULTS 
1.   Fidelity  Models 

The  available  signals  at  the  simulation  output  con- 
sist of  vectors  estimating  the  pilot's  perceived  motion  and 
attitude.   Each  human  dynamic  orientation  program  (see  Fig.  8) 
computes  perceived  linear  acceleration,  angular  velocity,  and 
attitude  vectors.   Algorithms   for  comparing  these  time  vary- 
ing signals  were  constructed  based  on  the  assumption  that  the 
greater  the  difference  between  corresponding  vectors  for  the 
aircraft  and  simulator,  the  less  the  fidelity.   The  following 
vector  variables  are  available  as  inputs  to  the  fidelity 
function.   (Notation  is  consistent  with  that  used  by  Young  in 
Ref.  12.) 

A^t)  =  perceived  aircraft  linear  acceleration 
Ag (t)  =  perceived  simulator  linear  acceleration 
w^t)  =  perceived  aircraft  angular  velocity 
Wgtt)  =  perceived  simulator  angular  velocity 
DOWN-,  (t)  =  perceived  aircraft  attitude 
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D0WN3(t)  =  perceived  simulator  attitude 

where:   a  =  aircraft,   s  =  simulator,  ^  indicates  a  perceived 

quantity 
(t)  indicates  a  time 
varying  function 

Three  different  fidelity  functions  were  constructed 

to  be  compared  with  the  experimental  data.   The  first  is  an 

average  percent  error  model  as  suggested  by  Young: 

"When  the  [motion]  sensations  are  clearly  supra- 
threshold,  the  most  likely  candidate  is  just  percent 
error,  the  ratio  of  perceptual  error  to  the  correct 
quantity.  When  the  model  indications  for  'correct' 
perceptions  are  subthreshold,  it  seems  more  logical 
to  assess  a  large  penalty  for  errors  that  are  large 
compared  to  the  threshold  value." 

For  each  sample  output  (one  per  second)  the  following  "cost 

to  realism"  indices  were  computed. 


A  A(t)| 

=  C, (t)  (2a) 


max 


f  i  ^     i  ^   *\     A 


|Aw(t)| 

T7 ^ — 'l=cw(t)  (2b) 

max(|wa(t)|/  wthrj 

|A  DOWN(t)[        =  CD(t)  (2c) 

where:   A  =  the  difference  between  aircraft  and 
simulator  vectors 

|x|  =  amplitude  of  vector  x 

thr  =  threshold  value 

Threshold  values  used  in  the  simulation  were  those 
computed  in  Ormsby's  vestibular  modelling  thesis  [Ref.  11] 
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Athr  =  -005  "G,s" 

w,,   =  .0146  radians/sec 

The  vector  DOWN  is  a  constant  amplitude  unit  vector;  thus, 
the  cost  index  is  simply  the  magnitude  of  the  difference 
vector.   Averages  for  each  of  the  three  cost  indices  were 
computed  over  each  simulator  run.   Though  the  task  was  de- 
signed to  last  three  and  one-half  minutes,  subjects  often 
used  over  four  minutes  to  complete  it.   Aircraft  acceler- 
ation data  from  each  run,  up  to  four  minutes,  was  used  as 
input  to  the  simulation.   The  following  average  cost  values 

were  computed: 

n 


i=l 


n 


^Cv^)  <3b> 


n 


cd  ■  "LcD(t  ) 


where:   n  =  number  of  output  time  frames 

Several  subjects  noted  that  during  many  runs  there 
existed  unexplained  and  unrealistic  motion  excursions  about 
the  rotational  axis,  but  more  pronounced  in  the  translational 
cues.   A  measure  using  the  same  cost  indices  which  is  more 
sensitive  to  these  rapid  anomalies  (if  caused  by  the  program- 
ming and  not  the  hardware)  is  a  root  mean  square  (RMS)  fig- 
ure.  The  second  fidelity  function  computes  RMS  cost  values 
as  follows: 
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c 


V  n   i=1  D   i 


(4a) 


(4b) 


(4c) 


It  is  not  unreasonable  to  assume  that  pilot  responses 
may  simply  be  proportional  to  the  amplitude  of  the  error  vec- 
tor.  The  third  fidelity  function  computes  the  average  of 
each  error  vector;  the  numerators  of  equations  3a,  3b,  and 
3c. 


n 


--1 


1 
w  =  - 

n 


D=  1 

n 


|AA(ti)| 


1=1 

n 


i=l 

n 


Aw(t.). 


A  DOWN(t. 


(5a) 


(5b) 


(5c) 


i=l 


A  total  of  133  simulator  runs  were  recorded  success- 
fully; 97  of  which  used  wash-out  parameter  sets  other  than 
the  standard  set.   Minimum,  maximum,  and  average  cost  indices 
computed  by  the  three  fidelity  functions  are  shown  in  Table 
VII. 
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MINIMUM 


AVERAGE 


MAXIMUM 


'w 


w 


W 


.9959 

1.0788 

1.4830 

.5776 

.7603 

.9457 

.0626 

.0898 

.1243 

.9983 

1.1428 

2.0106 

.6750 

.8309 

1.0671 

.0715 

.1045 

.1451 

.0403 

.0612 

.0911 

.0108 

.0179 

.0282 

.0626 

.0898 

.1243 

Table  VII:   Cost  Index  Values  Computed  by  the  Simulation 

2 .   Experimental  Optimization 

Response  of  the  simulation  to  changes  in  the  wash-out 
parameters  was  investigated  to  determine  its  adaptability  to 
experimental  optimization  methods  such  as  an  optimal  gradient 
search.   If  the  simulation  is  to  be  an  effective  research 
tool  for  determining  optimal  control  system  variables,  then 
its  response  characteristics  must  be  of  a  tractable  form. 
Time  did  not  allow  for  a  full  exploration  of  its  response 
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"surface"  by  an  optimization  routine,  nor  were  manual  methods 
employed  to  search  for  a  minimum  "cost"  solution. 

Linear  and  rotational  wash-out  filter  gains  were 
varied  over  the  range  of  the  parameters  used  in  the  experi- 
ment.  The  response  function  was  arbitrarily  chosen  to  be: 

f(CA,Cw,CD)   =   CA  4-  Cw  +  CD  (6) 

The  linear  gains,  Kt,  were  fixed  at  different  levels 
while  the  rotational  gains,  K  ,  were  varied  together  from 
0.0  to  0.8  at  intervals  of  .05.   The  response  curves  obtained 
are  shown  in  Fig.  10.   The  basic  parameter  set  used  during 
these  simulation  runs  was  Pq . 

C.   SIMULATION  VALIDATION 

Validation  of  the  simulation  consisted  of  comparing  the 
subjective  fidelity  ratings  with  predictions  from  the  simu- 
lation.  It  was  then  a  question  of  which  ratings  and  cost 
indices  to  analyze.   Since  the  strongest  correlations  be- 
tween subjects  were  for  ratings  of  the  combined  motion  cues, 
the  combined  ratings  were  compared  to  the  following  unweight- 
ed linear  combinations  of  the  cost  indices: 

Cfl  +  Cw  ♦  CD  (7a) 

CA  +  C   +  C^  (7b) 

A     w     D 

A   +  W   +  D  (7c) 

Spearman  rank  correlation  coefficients  were  computed  for 
each  subject  and  cost  index  combination.   (Fidelity  ratings 

56 


2.05  -- 


1.90  " 


0.0 


0.2  0.4  0.6 

Rotational  Gains,  K 


0.8 


Figure  10:   Simulation  Response  to  Varying  Filter  Gains 
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were  subtracted  from  a  constant  to  yield  scores  propor- 
tional to  the  cost  to  fidelity.)   The  results  are  shown  in 
Table  VIII. 


___^^^_     _     _N0. 

C,  +  C  ,  +  Cn    C,  +  C   +  C-    A   +   W   +   D    RUNS 
A     W      D      A     W      D 


S 
U 
B 
J 
E 
C 
T 
S 


.4147* 

.5138** 

.2673 

14 

.2893 

.3445* 

-.0403 

21 

.3802** 

.3986*** 

-.1489 

26 

.0542 

.2965 

-.0599 

13 

.3217* 

.3227* 

.3953** 

23 

Significance  levels:   *p<0.1,  **p<.05,  ***p<.025 


Table  VIII:   Spearman  Rank  Correlations  Between  Subjects 
and  Simulation  Cost  Indices 


Weighted  sums  of  the  cost  indices  were  computed  by  mul- 
tiple linear  regression  techniques  in  order  to  improve  the 
correlations.   Using  the  assumptions  for  parametric  methods 
given  for  the  analysis  of  variance  above,  a  regression  model 
of  the  form: 


RC  =  alCw  +  a2CA  +  a3CD  +  a4 


(3) 


where:   R^  = 


-   R, 


Rc  =  combined  motion  rating 
a   =  regression  coefficients 
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was  used  to  compute  regression  coefficients  for  each  subject 
and  for  all  subjects  combined.   Results  are  shown  in  Table  IX. 


Regression  Coefficients 


Sig.@ 
MULT.  R   Level 


S 

u 

B 
J 
E 
C 
T 
S 


5 

All 
Subjects 


-.530 

2.697* 

5.676 

-.874 

.5660 

.2571 

2.819 

1.845 

-17.628* 

.547 

.4926 

.1826 

3.578 

2.030 

1.811 

-2.926 

.3671 

.3541 

-6.540 

.449 

-7.916 

9.312 

.2949 

.8346 

8.792* 

.780 

9.096 

-4.999 

.4589 

.2030 

-1.336 

1.513** 

-3.070 

2.653 

.2442 

.1245 

@  Significance  level  associated  with  the  ANOVA  F  ratio  for  the 
regression  sum  of  squares 

Significance  levels:  *p<0.1,  **p<0.05 


Table  IX:   Regression  Coefficients  for  Combined  Ratings 
versus  RMS  Cost  Indices 
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IV.   DISCUSSION 

Results  of  the  experiment  tend  to  reject  the  hypothesis 
that  similar  motion  cues  were  presented  to  the  subjects.   As 
shown  in  Figure  9,  a  cluster  of  nine  parameter  sets  exists 
toward  the  center  of  the  range  of  means.   Parameter  sets  one 
and  three  appear  significantly  better  than  the  others,  while 
sets  9,  11,  14,  and  15  are  significantly  inferior  to  the 
others.   The  fact  that  the  center  group  of  parameter  sets 
produced  indistinguishable  motion  cues  is  a  probable  explana- 
tion for  the  weak  correlations  between  subjects.   Parametric 
Pearson  product-moment  correlations  were  computed  yielding 
similar  results  at  comparable  significance  levels.   The 
highly  significant  coefficients  of  concordance  also  indicate 
that  ordinal  agreement  between  subjects  was  less  than  antic- 
ipated.  The  overall  higher  correlations  for  the  combined 
motion  ratings  infer  that  subjects  were  capable  of  yielding 
more  reliable  scores  after  rating  both  linear  and  rotational 
motion  cues  on  an  absolute  scale,  followed  by  a  pairwise 
comparison  with  the  standard  motion.   The  higher  correlations 
could  also  be  explained  by  the  fact  that  several  subjects 
said  their  responses  resulted  from  an  overall  impression  of 
the  motion  cues  rather  than  a  subjective  average  of  the 
linear,  rotational,  and  comparative  scores.   Two  subjects 
rated  the  combined  motion  cues  higher  than  either  linear  or 
rotational  scores  for  several  runs. 
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Comments  by  the  pilots  indicate  that  a  large  percentage 
of  the  variance  of  the  pilot  ratings  is  likely  due  to  un- 
realistic dynamics  of  the  motion  drive  mechanism.   High 
frequency  vibrations  in  the  system  were  not  as  distracting 
and  did  not  seem  to  lessen  fidelity  as  much  as  the  low  fre- 
quency, larger  amplitude  anomalies  experienced  by  the  sub- 
jects.  Whether  these  were  caused  by  the  software  or 
hardware  is  not  known. 

Cost  indices  computed  by  the  simulation  show  large  dif- 
ferences between  outputs  of  the  vestibular  models.   Errors 
ranged  from  58  to  148  percent  of  the  perceived  aircraft 
motion  vectors  as  shown  in  Table  VII.   Whether  this  typifies 
the  capabilities  of  the  S.01  simulator  or  is  largely  a  result 
of  programming  limitations  and  inaccuracies  is  not  readily 
apparent.   A  significant  amount  of  error  is  most  likely  due 
to  roundoff,  precision  limits,  and  the  discrete  nature  of  the 
program.   Motion  histories  of  the  aircraft  and  simulator  have 
cycle  times  of  .05  seconds,  while  the  output  rate  of  each 
vestibular  model  is  one  per  second.   This  represents  a  con- 
siderable loss  of  information.   If  the  output  sampling  rate 
equaled  the  input  rate,  it  may  have  been  more  responsive  to 
the  rapid  motion  anomalies  noted  by  the  subjects.   A  quick 
excursion  lasting  less  than  a  second  may  not  have  been  re- 
flected in  the  output  of  the  human  dynamic  orientation  model. 

The  aircraft  model  used  during  the  experiment  was  that  of 
the  F-15  fighter.  Ideally,  subjects  with  recent  F-15  experi- 
ence should  have  been  chosen,  but  due  to  time  and  availability 
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constraints  selection  of  the  pilots  was  less  than  optimal. 
It  was  felt  that  this  factor  contributed  greatly  to  the  lack 
of  higher  correlations  of  the  ratings  between  subjects. 
This  conclusion  was  supported  by  the  weak  correlations  be- 
tween the  first  and  second  replications  of  the  parameter 
sets.   This  also  implies  an  inability  of  the  subjects  to 
repeat  the  task  as  well  as  inconsistent  responses  between 
subjects . 

Response  of  the  simulation  to  changes  in  wash-out  filter 
gains  indicated  a  definite  tendency  for  convexity  in  the 
neighborhood  of  minimum  output.   This  can  be  explained  in- 
tuitively.  When  the  gains  are  relatively  low,  the  lack  of 
motion  cues  produces  a  large  discrepancy  between  perceived 
aircraft  and  perceived  simulator  motions.   When  the  gains 
are   high,  cues  become  overly  rapid,  exagerated,  and  fre- 
quently cause   the  simulator  to  reach  its  "soft"  programming 
limits.   When  reached,  the  limits  force  the  linear  and  ro- 
tational cues  to  become  uncoordinated,  reducing  fidelity. 
The  simulation  output,  defined  by  equation  6,  varied  over  a 
small  range  as  the  rotational  gains  were  varied  as  shown  in 
Figure  10.   For  each  level  of  the  linear  gains,  the  output 
varied  only  2  to  4  percent.   Parameter  sets  which  produced 
significantly  better  motion  cues,  P1  and  P3,  were  comprised 
of  filter  gains  set  at  0.3.   The  tendency  for  the  simulation 
output  to  minimize  at  approximately  0.3  indicates  a  capacity 
to  predict  subjective  preferences  for  motion  cues  produced 
by  these  parameter  sets.   The  simulation  was  also  run  with 
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the  linear  gain,  K^,    set  at  0.   Response  was  less  than  that 
for  K1  set  to  0.1  throughout  the  range  of  K  .   This  implies 
that  optimal  motion  cues  include  very  little  linear  acceler- 
ation for  the  specific  flight  used  as  input  to  the  simulation. 

Comparison  between  the  simulation  output  and  subject 
responses  is  summarized  as  follows: 

■  Positive  correlations  were  obtained  for  the  first  two 
fidelity  models,  equations  7a  and  7b  (see  Table  VIII). 

■  The  third  function,  equation  7c,  showed  little  correla- 
tion with  the  exception  of  subject  five. 

■  The  second  fidelity  function  consistently  produced 
higher  and  more  significant  correlations  infering  that 
the  RMS  cost  indices  are  a  better  method  for  simulating 
subject  responses. 

■  Linear  regression  using  the  RMS  cost  indices  as  inde- 
pendent variables  and  the  combined  motion  cue  ratings 
as  dependent  variables (equation  8)  showed  little  con- 
sistency of  regression  coefficients  among  the  five 
subjects  (see  Table  IX) . 

■  Linear  regression  analysis  yielded  multiple  correla- 
tions not  appreciably  better  than  the  correlations 
for  the  unweighted  cost  index  combination. 

■  Very  little  correlation  was  found  between  ratings  of 
the  linear  and  rotational  motion  cues  and  the  corres- 
ponding cost  indices,  such  as  the  rotational  ratings 
and  Cw. 
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Nine  of  the  15  wash-out  filter  parameter  sets  were 
grouped  in  the  center  of  the  range  of  mean  responses  as 
shown  in  Figure  9.   Since  these  did  not  represent  signifi- 
cantly different  means,  Spearman  rank  correlations  were  com- 
puted between  the  subjects*  combined  scores  and  the  second 
fidelity  function,  equation  7b,  for  the  remaining  six  param- 
eter sets.   The  results  are  shown  in  Table  X. 


No.  of  Runs 


SUBJECTS 

sl 

S_2 

fa 

f£ 

fs 

710** 

.0748 

.6545* 

.5000 

.7660* 

6 

7 

10 

3 

8 

Significance  levels:  *p<.05,  **p<.01 


Table  X:   Spearman  Rank  Correlations  for  the  Remaining 
Six  Parameter  Sets 


With  the  exception  of  subject  two,  the  subjects  showed  a 
marked  increase  in  correlation  with  the  simulation  results. 
Although  the  sample  sizes  for  this  post  hoc  analysis  were 
small,  three  of  the  five  correlations  were  significant.  It 
was  noted  that  the  correlations  compared  closely,  rankwise, 
to  the  flight  experience  each  pilot  had  in  high  performance 
jet  aircraft  as  indicated  in  Table  I. 
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VI.   RECOMMENDATIONS 

Recommendations  for  further  research  include  suggestions 
for  improvement  of  the  experimental  data  base  and  the  com- 
puter simulation.   Better  reliability  and  consistency  in 
subject  responses  and  a  more  accurate  and  responsive  simu- 
lation would  improve  the  prediction  capability  of  the 
simulation. 

Further  studies  should  employ  subjects  proficient  in  the 
simulated  aircraft.   Since  the  results  appeared  to  be  strong- 
ly affected  by  the  subjects'  ability  to  repeat  the  task,  use 
of  experienced  test  pilots  familiar  with  flight  dynamics 
would  help  create  a  more  reliable  data  base.   Future  re- 
search should  include  preliminary  studies  to  determine  other 
control  parameters  which  when  varied  would  create  motion 
cues  more  significantly  different.   Discrete  rating  scale 
data  is  difficult  to  analyze  with  classic  analysis  of  vari- 
ance techniques,  therefore,  methods  for  obtaining  pilot  re- 
sponses from  a  continuous,  interval  scale  should  be  pursued. 

Time  did  not  allow  for  elimination  of  the  inefficiencies 
of  the  simulation.   Recomputation  of  the  human  dynamic  orien- 
tation program's  input  and  output  model  parameters  would  most 
likely  improve  the  responsiveness  of  the  model.   Converting 
the  program  to  double  precision  would  reduce  the  inaccuracies, 
but  double  the  computer  memory  requirements.   The  core  re- 
quired as  listed  in  Appendix  C  and  as  executed  on  the  IBM 
360/67  Operating  System  was  458K.   The  entire  flight 
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acceleration  data  set  was  read  in  at  once.   Other  time  his- 
tories, such  as  the  simulation  motion  vectors,  were  computed 
creating  large  core  requirements.   If  each  time  frame  were 
processed  through  the  simulation  before  the  next  was  read  in, 
the  core  requirements  would  be  reduced  by  7  5  percent. 

It  was  assumed  that  an  accurate  facsimile  of  the  simu- 
lator motion  was  generated  by  the  drive  program.   Since  the 
lead  compensation  network  was  bypassed  in  the  simulation, 
follow-on  studies  should  compare  the  output  of  the  drive 
program  with  actual  simulator  motion  histories.   This  is  a 
possible  source  of  considerable  error. 

The  motion  vectors  for  the  simulator  and  aircraft  could 
have  been  compared  directly  without  use  of  the  human  dynamic 
orientation  program.   A  similar  analysis  of  the  difference 
between  these  vectors  would  provide  a  way  of  validating 
Ormby ' s  model.   If  the  method  incorporating  the  model  cor- 
related appreciably  higher  with  subject  ratings,  a  reasonable 
validation  of  the  human  dynamic  orientation  would  be  achieved. 
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APPENDIX  A 
INSTRUCTIONS  TO  SUBJECTS 

A.  PURPOSE  OF  THE  EXPERIMENTS 

The  intent  of  this  study  is  to  analyze  pilot  response  to 
varying  degrees  of  simulator  realism  or  fidelity.   Subjective 
ratings  by  commercial  pilots  like  yourself  will  help  con- 
struct a  mathematical  model  of  the  human  process  of  compar- 
ing simulator  motion  characteristics  with  those  of  the 
simulated  aircraft.   Ultimately,  our  goals  are  to  validate  a 
model  of  human  dynamic  orientation  and  to  optimize  the  drive 
programs  of  aircraft  simulators. 

B.  TASKS 

You  will  be  asked  to  repeat  two  familiar  tasks,  each 
lasting  from  three  to  three  and  one  half  minutes.   The  first 
will  be  a  basic  airwork  pattern,  designed  to  employ  the 
maneuvers  common  to  operations  in  the  terminal  phase  of 
flight.   The  second  task  is  a  full  ILS  approach  to  a  wave- 
off.   The  basic  airwork  task  will  be  repeated  approximately 
fifty  times,  while  the  drive  program  parameters  are  varied 
to  produce  different  degrees  of  fidelity.   At  the  end  of 
each  task  you  will  be  asked  to  rate  the  characteristics  of 
the  motion  using  three  rating  scales  discussed  below,  as 
well  as  provide  your  subjective  opinion  of  the  motion,  task, 
etc.   During  subsequent  sessions,  you  will  rate  the  motion 
following  the  ILS  approach  task,  repeated  roughly  the  same 
number  of  times. 
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Each  session  will  consist  of  approximately  12-15  flights 
and  last  no  longer  than  80  minutes.   Detailed  diagrams  of 
each  task  will  be  available  during  the  runs.   Additional 
details  will  be  provided  before  the  sessions  begin. 

The  simulated  aircraft  in  this  study  is  the  F-15,  a  high 
performance,  fighter  aircraft.   Although  you  probably  have 
had  no  experience  in  the  F-15  itself,  your  background  in 
high  performance  type  aircraft  was  a  major  consideration  in 
your  being  selected  as  a  pilot  for  this  study. 

C.   USE  OF  THE  RATING  SCALES  -  OBJECTIVE  RESPONSES 

Each  trial  run  will  employ  a  different  set  of  drive  pro- 
gram parameters.   A  standard  set  of  parameters  will  produce 
a  "standard  motion"  which  will  be  presented  during  the  first 
run  and  every  fourth  run  thereafter.   Repetition  of  the 
standard  motion  will  reinforce  your  perception  of  it,  allow- 
ing you  to  make  comparative  evaluations  of  the  trial  runs. 
A  typical  sequence  of  12  runs  per  session  will  be  presented 
as  follows: 


Trial  No. 


7    8 


10    11    12 


where   5 


indicates  that  the  fifth  run  of  the  session  will 
be  the  "standard  motion." 

1.   Component  Rating  Scale 

Evaluate  the  linear  (translational)  and  rotational 
components  of  the  motion  following  each  simulator  run  using 
the  component  rating  scale.   Rate  the  motion  characteristics 
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based  on  the  descriptive  phrases  provided  at  each  level  of 
the  scale.   These  subjective  descriptions  should  be  used 
merely  as  an  aid  in  your  evaluations.   Responses  should 
reflect  consideration  of: 

a.  Response  to  Control  Inputs  -  simulator 
reaction  to  movement  or  pressure  on  flight  controls. 

b.  Motion  Deficiencies  -  motion  deviations  from 
experienced  or  anticipated  aircraft  motion,  not  a  result  of 
control  inputs  (e.g.  response  to  turbulence,  oscillations, 
etc.  )  . 

c.  Other  -  any  other  discrepancies  in  the  motion 
dynamics  which  you  detect  and  consider  significant. 

Your  responses  should  be  formulated  in  reference  to 
the  absolute  scales  as  well  as  in  comparison  to  your  rating 
of  the  standard  motion. 

2 .  Comparative  Rating  Scale 

Rate  the  trial  runs  in  comparison  the  the  standard 
motion  using  the  comparative  rating  scale.  Brief  descrip- 
tions of  each  level  of  the  scale  should  aid  your  evaluation. 

3.  Combined  Rating  Scale 

Using  the  descriptive  phrases  on  the  component  rating 
scale,  rate  the  trial  run  on  an  overall  or  combined  basis. 
Reference  should  be  made  to  your  scores  given  on  both  the 
component  and  comparative  rating  scales. 

4 .  General 

a.   If  you  cannot  choose  between  two  adjacent  levels 
on  either  of  the  scales,  a  "half"  level  score  is  allowed, 
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provided  careful  consideration  is  given  to  both  adjacent 
levels. 

b.  Ensure  that  your  rating  of  the  "standard  motion" 
is  recorded  on  the  scale  and  referred  to  while  rating  trial 
runs. 

c.  Basic  performance  measures  will  be  recorded 
only  during  the  ILS  approaches.   Although  your  performance 
will  be  compared  to  that  of  other  pilots,  our  interest  lies 
in  how  your  performance  varies  with  changes  in  fidelity. 

d.  Try  to  ignore  changes  in  your  performance  on 
the  task  as  each  session  progresses. 

e.  Your  evaluations  should  reflect  your  assessment 
of  all  the  motion  cues  detected. 

D.   SUBJECTIVE  COMMENTS 

After  rating  each  trial  run,  brief  comments  on  the  items 
listed  below  are  encouraged.   It  is  not  necessary  to  remark 
on  each  issue,  but  only  on  those  you  feel  significant  or 
particularly  relevent  to  the  trial  run. 

1.  Your  ability  to  rate,  separately,  the  linear  and 
rotational  components  of  the  motion. 

2.  Subjective  description  of  the  motion  characteristics. 

3.  Difficulties  encountered  with  the  task  as  it  is 
affected  by  the  particular  motion  characteristics. 

4.  Impressions  of  the  evaluation  task. 

5.  Duration  of  the  task  versus  your  ability  to  evaluate 
the  motion. 
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6.  Physical  discomfort  (fatigue,  nausea,  etc.). 

7.  Workload  level  in  terms  of  the  attention  level 
required  to  perform  the  task. 
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APPENDIX  B 

WASH-OUT  FILTER  PARAMETER  SETS 

Each  wash-out  filter  (one  for  each  degree  of  freedom) 
is  described  by  three  parameters.   The  values  of  the  param- 
eters used  in  the  experiment  are  shown  in  Table  Bl.   Where 
only  one  figure  is  shown,  that  value  was  used  for  all  three 
channels.   Natural  frequencies  for  the  translational  chan- 
nels are  in  the  order  of  roll,  pitch,  and  yaw.   The 
parameters  are  defined  as  follows: 


K   -  Linear  gains  for  the  translational  channels 

Kr  -  Linear  gains  for  the  rotational  channels 

w   -  Natural  frequencies  for  the  translational  channels 

w   -  Natural  frequencies  for  the  rotational  channels 

r 

£,  -  Damping  factors  for  the  translational  channels 
£   -  Damping  factors  for  the  rotational  channels 


The  "standard  motion"  parameter  set  is  denoted  by  P 


0 


72 


Parameter 

Set  No.     Kn       K  w,        wv 

1       r  1        r 


p10 


P12 


P13 


P14 


P15 


k  £r 


PQ  0.5  0.5  .4  .4  1.0  1.2  .707  .707 

Px  0.3  0.3  .2  .2   .2  0.7  .707  .707 

P2  0.1  0.1  .4  .4  1.0  1.2  .707  .707 

P  0.3  0.3  .2  .2   .7  1.6  .707  .707 

P4  1.0  0.1  .7  .7  1.5  0.7  .707  .707 

P5  1.0  0.1  .7  .7  1.5  1.2  .707  .707 

P^  1.0  0.1  .7  .7  1.5  1.6  .250  .707 
6 

P  0.7  0.3  .7.71.5  1.2  .707  .707 

P  0.1  1.0  .7  .7  1.5  1.6  .707  .707 
8 

P  0.3  0.7  .2  .2   .7  0.7  .707  .707 


0.3      0.7      .2  .2   .7     1.2      .707     .707 


p  0.1      1.0      .7  .7  1.5     0.7      .707     .707 


0.1  1.0  .7  .7  1.5  1.2  .707  1.5 

0.3  0.7  .7  .7  1.5  1.6  .707  .707 

0.7  0.7  .7  .7  1.5  0.7  .707  .707 

1.0  1.0  .4  .4  1.0  1.2  .707  .707 


Table  Bl:   Wash-out  Parameter  Sets  Used  in  the  Experiment 
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APPENDIX  C 
SIMULATION  SOURCE  PROGRAM  LISTING 

The  simulation,  written  in  FORTRAN  IV  H,  consists  of  sub- 
routines which  perform  the  functions  shown  in  the  block  dia- 
gram in  Figure  9.   The  program  listing  in  this  appendix  is  in 
the  configuration  used  to  perform  the  experimental  optimiza- 
tion described  in  section  III.B.2.   The  main  subroutine 
(called  SIM)  controls  the  other  subroutines  and  prints  the 
wash-out  parameter  set  and  messages  displaying  the  progress 
of  the  program  execution.   A  main  program  (not  listed)  reads 
in  the  values  for  the  wash-out  parameters,  varies  them,  and 
iteratively  calls  the  main  subroutine.   The  flight  acceler- 
ation time  histories  are  read  in  by  subroutines  READ.   The 
model  parameters  used  by  the  human  dynamic  orientation  pro- 
gram (subroutine  HDO)  are  those  obtained  with  the  program 
listing  contained  in  Ref.  21.   The  simulation  was  executed 
on  an  IBM  360/67  Operating  System  using  single  precision, 
floating  point  variables. 

Subroutines  MCPY  and  MINV  are  IBM  system  library  sub- 
routines called  from  subroutine  COTRN  to  perform  matrix 
operations.   MCPY  copies  one  array  variable  into  another. 
MINV  inverts  a  square  matrix. 
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*  MAIN  SUBROUTINE  * 


C 

c 
c 
c 


c 
c 
c 

IOC 

c 


PROGRAM  TO  EVALUATE  THE  SOFTWARE  PARAMETERS  OF  THE  S.Ol 
*CTICN  SIMULATOR  USING  THE  ORMSBY  HUMAN  DYNAMIC 
ORIENTATION  MODEL 

SLBROUTINE  SI M ( CAD, CA VG,CRMS , K ) 

COMMON/ XFLOAT/TMT (3,3) , TMR ( 3,3) ,BTACC(3 ) , BRACC(3) 
*,  TDDS(3) ,RSFX(3 ) 

CCMMON/IFIXED/IMODE, IDUMFU, I  FORD,  ICHG 

CCMM0N/WSH0UT/AKR(3)  ,AKK  3) ,OMGR( 3),OMGT (3  )  ,ZETR(3) 
*,  ZETT(3) 

CCMMCN/FLIGHT/XA1 4300,3) , XR ( 48  00 , 3) , P, C , IP AR 

DIMENSION  ACC(3 ) ,R0T(3 ),ANG(480  0, 3) 
DIMENSION  SFS<4800,3) ,WS( 48  00, 3) , WA(4800,3) 

CATA  DT/. 05/, N/4800/ 

PRINT  THE  WASHOUT  FILTER  PARAMETERS 

WRITE(6,100)  _   „ 

FORMAT',  1H0,  'WASHOUT  PARAMETERS  READ  FROM  DATA  FILE:*) 

WRITE(6,120  )(AKR(  I ) , I=L ,3) 


C 
C 

C 


v*RITE(  6  ,121  MAK1  (  I) 
WRITE(6, 122)(OMGR(I 
/sRITE(6,123)(0MGT(  I 
*RITE(6,124)(  ZETRd 
WRITE(  6,  125H  ZETTd 


1=1,3) 

,1=1,3) 

,1=1,3) 

,1=1,3) 

,1=1,3) 


REAC    THE    FLIGHT    DATA     (AIRCRAFT    MODEL    OUTPUT) 

IF(K.GT.l  )    GO    TO    131 

CALL    REAO(NSP) 
C 

109    FCRMATt  •   FLI3HT  OATA  READ  FROM  TAPE  FILE') 

r 

C        INITIALIZE    THE    INPUT    ACCELERATIONS    TO    ZERO 

C 

131  BT£CC(1  )=0. 

BTACC(2)=0. 

ETACC( 3)=-32.2 

ERACC( 1 )=0. 

ERACC(2)=0. 

ERACC(3 )=0. 


C 
C 
C 


101 

c 
c 
c 


SET     THE    "FLAGS"    FOR     SUBROUTINE    S01MGT 

IQUMRJ=0 

IFCRD=0 

ICFG=1 

SET    THE    MODE    TO    "INITIALIZE"    AND    CALL    SUBROUTINE    S01MCT 

IMCDE=-1 

CALL    SOIMOT(DT) 

FCRMATP  ^SUBROUTINE    S01M0T     INITIALIZED') 
SET    THE    MODE    TO    "RUN" 

IMCDE=1 
SET    THE     INPUT    ACCELERATION    SAMPLE    AND    CALL     SUBROUTINE 
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C       S01MCT    TO    CREATE    THE    DRIVE    COMMANDS 

DC    20    1=1, NSP 
C 

ETACC( 1)=XA(I, 1) 

BTACC(2)=XA(I,2) 

8TACC(3)=XA(I,3) 

2RACC( 1 )=XR(I,  I  ) 

ERACC(2 )=XR(I,2  ) 

8RACC(3)=XR(I,3) 
C 

CALL    SOIMOT(DT) 
C 

C       TRANSFORM    THE    DRIVE     SIGNALS    BACK    INTO    CA<3    (PILCT)    COORDI- 
C       NATES 
C 

CALL    C0TRN(TDDS,TMT,1,ACC ) 

CALL    COTRN( RSFX,TMR,-1 ,ROT) 
C 

SFSd  ,1  )=ACC(1) 

SFSd  ,2)=ACC<2) 

SFS( 1,3 )=ACC(3)-32.2 

ANG( I ,1 )=ROT( 1 ) 

ANG< I, 2 )  =  ROT(  I) 

ANG( I,3)=RCT(3) 
C 

20  CONTINUE 

C 

V»RITE<6,102) 

102  FCRMAT(«        SIMULATOR    DRIVE    SIGNALS    CREATED') 
C 

C       DIFFERENTIATE    THE    ROTATIONAL    POSITION    COMMANDS    TO    08TAIN 
C       HCTATICNAL    VELOCITY     FOR    USE    BY    THE    HUMAN    DYNAMIC 
C       ORIENTATION     (HOC)     SUBROUTINE 
C 

CALL    DIFF<DT,ANGf WS,NSP»NJ 
C 

WRITE(6,103 ) 

103  FCRMAT(  •   ROTATIONAL  DRIVE  SIGNALS  DIFFERENTIATED') 
C 

C   CALL  SUBROUTINE  HDO  TO  COMPUTE  THE  PERCEIVED  SIMULATOR 

C   MOTICN  VECTORS 

C 

CALL  HDC(DT,NSP,SFS,WS,9,N) 
C 

WRITE(6, 104) 

104  FCRMAT{«   PERCEIVED  SIMULATOR  MOTION  VECTORS  CREATED1) 
C 

C   INTEGRATE  THE  ROTATIONAL  FLIGHT  ACCELERATION  TO  BE  USED 

C   BY  SL3P0UTINE  HDO 

C 

IFtK.GT.l )  GO  TO  132 

CALL  INT(DT,XR,WA,NSP,N) 
C 

wRITE(6 ,105 ) 

105  FORMAT!'   ROTATIONAL  FLIGHT  ACCELERATION  INTEGRATED') 
C 

C   CALL  SUBROUTINE  HDO  TO  COMPUTE  THE  PERCEIVED  AIRCRAFT 

C   MCTICN  VECTORS 

C 

CALL  HDC(DT,NSP,XA,WA,8,N) 
C 

WRITE(6,107 ) 
107         FQRMAT(»       PERCEIVED    FLIGHT    MOTION    VECTORS    CREATED') 
C 

C       CPERATE    THE    COST    ANALYSIS    ROUTINE 
C 

132  CALL    COST(CAD,CAVG,CRMS) 

C 

wRITE(6,108) 
10*    FCRMAT(1H0,'  COST  ROUTINE  COMPLETE') 
C 
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c 
c 

120 
121 
122 
123 

124 
125 
C 


FCRMAT    STATEMENTS 


FGRMAT( /5X»  'AKR( 
FCRMAT(5X,»  AKK  1 
FCRMAT( 5Xt 
FCRMAT <  5X, 
FCPMAT( 5X, 
F0RMAT( 5X, 


GMGR( 

0MGT( 


1,2,3 
t  2,3) 


RETURN 
ENO 


—    i 


1  ,  <■ ,  o  ;  = 

1,2,3)  = 

ZETR(1,2,3)  = 

ZETT(lr2,3)  = 


,3F6.3) 

3F6.3) 

3F6.3) 

3F6.3) 

3F6.3) 

3F6. 3/) 


* 


-■-   J,      A,    -*-     ^.wWJ^sC^vi.-^^.^J,     J,   ^  J,    .t    X   X  Jr   vU  ■*.    J,   ^    O-    *l,   J,     -V     -  *-   -"*  •"*-    ^-   -''    ^»-    -''   -»-'   *-V    ~~    -'- 

-i*  -t*  -v  ^*  TT-irT*  *r  •t*  **»  ^r  -in  ir  «**  t*  *n  ^  r  ttt'^T'^'T  t*  tt  ^r  -nr*  ^n  n*  -,i*  ***  *c  -»*  -t*  -^  -^  *t» 

Op 

SUBROUTINES  * 


^^^^^^^^^^^^t^^^^^^^^^^^^^^^^^^^^^*^^^*^^^^^^^^^^^^^^^^5!1^15!15^^ 


SLBRCUTINE    TO     INTEGRATE    THE    TIME    SERIES    CONSISTING    OF    THE 
RCTATIONAL     FLIGHT    ACCELERATION     (AIRCRAFT    MODEL    OUTPUT) 

SUBROUTINE     INT  {  H  ,  Y  ,  Z  ,NSP  ,N  ) 

DIMENSION    Y(Nt 3) ,Z(N,3) 

HT=.3333333*H 
DC    7    J =1,3 
SLM1  =  Y( 2,J)+Y(2, J  ) 
SLM1=SUM1+SUM1 

SUM1  =  HT*(  Y(  1,  J)+S'JM1  +  Y(  3,  J)  ) 
ALX1=Y( 4,J)+Y(4,  J) 
ALX1=AUX1+AUX1 
ALX1=SUM1+HT*(Y(3 
AUX2=HT*( Y(  l,J)+3 
*Y(4tJ))+Y(6fJ) ) 
SUM2=Y( 5, J )+Y(5,J  ) 
SUM2=SUM2+SUM2 
SLM2=AUX2-HT*( Y(  4- 
Z(1,J)=0. 
AUX=Y(3 ,J) +Y(3,J) 
AIX=AUX+AUX 
Z( 2,J)=SUM2-HT-(Y  (2 
Z( 3,J)=SUM1 
Z(4,J)=SUM2 


J)+AUX1+Y( 5, J)  ) 
875*(Y(2? J)+V (5,  J)  )+2 


J)+SUM2+Y(6,J)  ) 


JJ+AUX  +  YU,  J)  ) 


625*(Y(3, J  )  + 


DO    A    I=7,NSP,2 

SLM1=AUX1 

SUM2=AUX2 

ALX1=Y<  1-1, J)+Y(  1-1  ,  J) 

ALX1=AUX1+AUX1 

AL'X1=SUM1+HT*(  Y(  I-2,J)+AUX1+Y(  I  ,J)  ) 

Z( 1-2, J )=SUM1 

IF( I-NSP) 3,6,6 

AUX2=Y( I ,J )+Y( I, J) 

AUX2=AUX2+AUX2 

AIX2=SUM2+HT*(Y( I - I , J ) + AUX 2+Y ( 1+1, J) ) 

Z( 1-1, J )=SUM2 

CONTINUE 

Z<NSP-1 t J)=AUX1 

Z(NSP, J)=AUX2 

GC    TO    7 

Z(NSP-1 ,J) =SUM2 

Z(NSP, J )=AUX1 

CONTINUE 

RETURN 
END 
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c 
c 

c 

c 
c 
c 
c 
c 
c 
c 


SUBROUTINE 
TFE  FLIGHT 


TO  PRODUCE  THE  SIMULATOR  DRIVE  SIGNALS  FROM 
DATA  OBTAINED  FROM  THE  MAIN  PROGRAM 


C 
C 

10 

c 


SLBROUTINE    SOIMOT(DT) 

MOTION    DRIVE    SYSTEM    FOR    S.01     (SIX-DEGREE) 


PREPARED    BY    J.DOUVILLIER    (NASA) 
PROGRAMMED    BY    D.ASTILL  (CSC) 

MODIFIED    BY    DATA    ANALYSIS    GROUP 


(CSC) 


3 
12 


SEPT     1975 

SEPT     1975 

DEC     1976 


CCMM0N/XFL0AT/TMTO,3) »TMR( 3,3) ,BTACC(3) ,BRACCO) 
:,     TDDS(3) ,RSFX(3 ) 

CCMMON/IFIXED/IMODE,  IDUMF'Jt  IFORD,  ICHG 

COMMON/ WSHOUT/AKR (3) , AK 1( 3) , OMGR ( 3 ) ,OMGT(3) ,  ZETR( 3) 

»     ZETT(3) 

3)  ,OMGRQ(3)  .CMGREO)  .OMGTEO) 
ZDLF(  3)  ,GNRS  (3),AKRL(3)  ,SO,MGR(  3) 
tTOMGRE (3) ,  0MGRE2O) ,0MGT2O) 
,  S0MGT2O) ,FOMGT3(3) ,0MGTE4(3) 


CCMMON/MOTION/OMGRL( 
It  ZETAQO) , ZETRLO), 
2,  ZOMGR (3),OMGQSQ(3) 
3t  TOMZK3)  ,FCMGTE(3) 
<*,     CMGTQO)  ,GLEAD(3  ) 


, FNCHTPO) 

,TDCSP(3) 

KFLAGO) 


DIMENSION    RPOSL  (3  ),TP0SL(3)  ,TMS<  3),RAC<  3)  ,TACO) 
1  ,FNOT(3),TDCSI( 3)  ,TDCS I I( 3 ) , FNCHI (3) , FNCHI I  (3 ) 
2,RHFI ( 3)  ,RHFI  H  3) , RDHI ( 3)  t RCH ( 3 )  ,RDH( 3 )  t FN A ( 3 ) 
3,RHFF3(3),ROLIO),RLF(3),RLO),RSFO),RDM(3  ),ROCH{  3) 
^►,FCF(3)  rRHF(3)  ,RS  (3)  .RICO)  ,TICO),RSD(3  ),FNCLO) 
5, TDSC(3 I  ,FCHI<  3)  ,FKN(  3)  ,FNASO)  ,FNCH(3) 
6,TDDH(3 ),TDDT(3 ) , ADOL ( 3 ) , TDS ( 3 ) ,TSIL( 3) 
7,TDDA(3) ,TDDG(3) ,TSD(3i ,AFU(6) , ADSD(6) 
8,FOLUPS (6),GNH(3),AIC(6) 

EQUIVALENCE    ( FCLUPS ( 1 ) , AF J ( 1 ) ) 
EQUIVALENCE    ( ADSO < 1 ) , TSD( 1 ) ) 

DATA    G/32.2/,RG/.  031056/ 

CATA    TDDG/0., 0. ,-32.2/ 

DATA    FNOT/3*!./ 

CATA    TIC/3*0./,RIC/3*0./ 

DATA    AJDL/5.6 ,6 .8,5.5/ 

DATA    GNH/3M.0/ 

DATA    AKXL/1./ 

DATA    AKYL/1 ./ 

DATA  RP0SL/3*.5326/ 

CATA  TP0SL/3*8./ 

DATA  NORES/0/ 


XCOS(X) =1 .-<X*X)/2. 
XSIN( X)=X-( X**3) /6. 

IF(  IMODE) 10,60,100 


DC  20    IN=1,3 

KFLAG( IN)=0 
^GLUPS( IN)=0. 
FGLUPS (  IN  +  3  )  =  0, 
ACSD(  IN  )=0. 
ACSD( IN+3)=C. 
RLF( IN)=0. 
TSC(IN) =0. 
TSIL( IN)=0. 
RCH( IN)-0. 
RL(  IN)=0. 
RCHK  IN  )=0. 
RHFF6( IN)=0. 
DTP=0. 

CPHSF=1 . 
CTHSF=1 . 
CPSSF=1 . 


FAST  SINE  AND  COSINE  APPROX 


CLEAR  LIMIT  FLAGS,  INITIALIZE 
FILTERS,  AND  ZERO  CRIVES 


78 


c 

c 

c 
c 


2C 


24 
25 


C 
100 


SPhSF=0. 
STHSF=0. 
SPSSF=0. 
CPHLF=1 . 
CTHLF=1 . 
CPSLF=1 . 
SPHLF=0. 
STHLF=0. 
SPSLF=0. 

TMS(IN)=.5/ADDL(  IN) 
CCNTINUE 


IF(  ICHG.EQ.l)    GO    TC    25 

IF(DTP-DT)     24,100,24 

ICHG=1 

CTP=OT 

CALL    SOIMIC(DTP) 

ICHG=0 


UPDATE  FILTERS  IF  ICHG=1 
OR  DT  CHANGED 


C 

C 
C 

c 
c 


c 
c 
c 


c 
c 
c 


51 

53 


2CC 


DC  53  J=l,6 

IF(  IDUMFU.EQ.O)  GO  TO  51 

AFIM  J)  =  ADSD(J) 

FCL'JPS(  J)=AOSD(J) 

GO  TO  53 

AFU( J)=FOLUPS(J) 

CCNTINUE 

FNATS=0. 

OC  200  IR=1,3 

FNA(IR)=  3TACC( IR )*RG 
FNAS( IR )=FNA( IR)*FNA( IR) 
FKN(  IR)=AK1  (IR) *FNA< IR) 
FNATS=FNATS+FNAS ( IR) 

RDM( IR) =AKP (IR)-BRACC ( IR) 
RDCH(  IR  )=RDM(  IR)-ZOMGR( IR) 

RCF(IR)=RCH(IR)+RDCH(  IR)*QT 

CCNTINUE 


COMPUTE  FOLLOW  UP  IF  NO  CA3  MOTION 


COMPUTE  LCAD  FACTORS 

ROTATIONAL  WASHOUT 
KCH(  IR)-SOMGR( IR  )*RHFFB( IR) 


RGTATIONAL  (CAB  TO  GIMBAL) 
TRANSFORMATION 
MATRIX 


TMR( 1,1  )=1, 

TMR( 1 ,2 )=-CPHSP*SPSSF/CPSSF 

TMR(1,3)=SPHSF*SPSSF/CPSSF 

TMR(2,1 )=0. 

TMR( 2,2)=CPHSF/CPSSF 

TMR(2,3 )=-SPhSF/CPSSF 

TMR(3,1)=0. 

TMR(3,2 )=SPHSF 

TMR(3,3 )=CPHSF 

RCF( 1)  =  RCH(  1  )+RCH( 2)*TMR(  1, 
RDH( 2)=RCH( 2)*TMR(2,2 )+RCH( 
RDH( 3)=RCH( 2)*TMR( 3,2)+RCH( 

FNAT=SQRT(FNATS) 

IF(FNAT.EQ.O.)     FNAT=1.0 

FNMXL^FKNU  )*AKXL/FNAT 
FNMYL=FKN(2 )*AKYL/FNAT 
FNMZL       =-SQRT( 1 . - FNMXL*FNMXL-FNMYL*FNMYL ) 

LCW    FREQ    RESIDUAL    TILT    RATE 
RDLK1  )=<-FNMXL*<  SPSL  F*CTriL  F  )/  C  PSLF    -FNM  YL  *CPHl_F*CTHLF 


2)+RCH( 3)*TMR( 1,3) 

3)*TMR(2,3) 

3)*TMR(3,3) 
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c 

c 


c 
c 
c 


c 
c 


*  +FNMZL*SPHLF*CTHLF) 

POL  I  (2)=(FNMXL*CTHLF/CPSLF    +  FNMY  L*SPHLF*STHLF 

*  +FNMZL*CPHLF*STHLF) 
RDLI(3)=(-FNMXL*SPSLF*STHLF-FNMYL*CPHLF*CPSLF*STH|_F 

*  +FNMZL*SPHLF*CPSLF*STHLF) 

DC  250   IR=1,3 

RHFFB( IR)=RHFF6( IR)+ROH( IR)*DT 

IF( IMODE.LT.O)  ROLI(IR)=0. 

RL(IR)  =  RL(IR) +RDLI  (IR)*OT 

RDLR=ZDLF( IR) *( RL ( IR ) *AKRL ( IR)-RLF(IR)) 

RLF(IR)=RLF(IR)+RDLR 

IF(NORES.EQ.l )RLF ( IR)=0. 

RDHH  IR)  =  RDHI(  IR)+RDH(  IR)*DT 

IF(IMOOE.GE.O)  GO  TO  245 
RDHK  IR  )=RDH(  IR) 
RHFI(IR)=0. 

RFFIK  IR)=( ROHI ( IR )*GNRS( IR ) ) /0MGRE2 ( I R ) 
245  CONTINUE 

RHF( IR  )=RDHI(  IR)*GNRS( IR)-TGMGRE( IR)*RHFI( IR) 

*  -0MGRE2(  IR)*RHFI I (  IR ) 
RHFK  IR)  =  RHFI(  IR)+RHF<  IR)*DT 
RFFIK  IR)=RhFII(  IR  )+RHFI(  IR  )*DT 

ROTATIONAL  POSITION  LIMIT 

RSF(IR)=RLF(IR)+RHF(IR)+RIC( IR) 
RSFX( I R ) =R  S F ( IR ) 

I F(ABS( RSFX (IR) ) .GT.RPCSL( IR) )  RSFX( IR  )  =SI  GM  RPOS  L  ( IR) 
£  RSFX(IR)) 

ROTATIONAL  POSITION  COMMAND  OUTPUT 

RSD( IR) =RSF(IR)+ROH( I R) /OMGRQ( IR) 

IF(A6S  <RSD<  IR)  )  .GT.RPGSLI  IR  )  )     RSDUR)=SIGN(  RPOSL<  IR) 
*,RSD( IR)  ) 


250    CONTINUE 


C 
C 
C 


CPHLF= 
SPHLF= 
CTHLF  = 
STHLF= 
CPSLF= 
SPSLF= 
CPHSF= 
SPHSF= 
CTHSF= 
SThSF= 
CPSSF= 
SPSSF= 


TMT( 1, 
TMT( 2, 
TMT(3, 
TMT( 1, 
TMT(2, 
TMT( 3, 
TMT(  I, 
TMT(2, 
TMT(3T 


XCOS( 
XSIN( 
XCCS( 
XSIN( 
XCOS( 
XSIN( 
XCOS( 
XSIN( 
XCOS( 
X  S  I N  ( 
XCOS( 
XSIN< 


RLF(  I 
RLFd 
RLF(2 

RLF(2 
RLF(  3 
RLF(  3 
RSF(  1 
RSF(  1 
RSF(  2 
RSF(2 
RSF(3 
RSF(3 


(3) 


1  )=CP 
1  )  =  SP 
)=-C 
)=-C 
)=CP 
)=CP 
)=SP 
)  =  -S 
)=-S 


SSF*C 

SSF 

PSSF* 

PHSF* 

HSF*C 

HSF*S 

HSF*S 

PHSF* 

PHSF* 


THS 


TRANSLATICNAL     (CAS    TO    LINEAR 
TRANSFORMATION 
MATRIX 
F 


DRIVE) 


STHSF 

SPSSF*CTHSF+SPHSF*STHSF 

PSSF 

PSSF*STHSF+SPHSF*CTHSF 

PSSF*CTHSF+CPHSP*STHSF 

CPSSF 

SPSSF*STHSF+CPHSF*CTHSF 


FNCL( 1)=STHLF*CPSLF 

FNCL( 2)=-CPHLF*STHLF^SPSLF-SPHLF*CTHLF 

F.\CL(3)=SPHLF*STHLF*SPSLF-CPHLF*CTHLF 

<=CHI(  1)=FKN(1 )-FNCL(l ) /AKXL 
FCHK  2)=FKM(2)-FNCL(2)/Ai<YL 
FCHK  3)  =FKN(3) 
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c 
c 


c 
c 
c 
c 
c 


: 


c 

c 
c 
c 


OC  270  IT  =  1  ,3 

IF( IMODE.GE.O)  GO  TO  267 
FNCHI(IT)=0. 

FNCHIK  IT)=(FCHI(  IT)*GNH(  IT) )/CMGT2< IT) 
267  CCNTINUE 

FNCH(IT)=GNH{ IT ) *FCHI ( I T) -TOMZT { IT)*FNCHI< IT)- 
*  CMGT2( IT)*FNCHII( IT) 
FNCHK  IT)=FNCHI(  IT)+FNCH(  IT)*DT 
FNCHIK  IT)=FNCHI  I{  IT)+FNCHI(  IT)*DT 
TDOH( IT)=G*(FNCH< IT)+FNCL< IT)) 

270  CONTINUE 

DO  300  I T= 1 , 3 

TCDA( IT)=0. 
DO  280  JT=l,3 


TCDA( IT )=TDOA( IT ) +TDDH( JT) *TMT( IT,JT) 
230  CONTINUE 

TDDT( IT)=(TDDA( IT)-TDDG( IT) )*CNOT(IT) 


DUMMY    FOLLOW-UPS 


IF( IDUMFU.EQ.L)     F0LUPS1 IT)=TSD( IT) 
C 
C  CHECK    OPERATING    AREA PARA3CLIC    LIMIT? 


IF(FOLUPS( IT)*TSIL(  IT) )     233,283,28  2 
2  92     ZTEMP  =  TSIL(  IT ) *TS I L(  IT ) *TMS I  IT ) 

IF(  (TPOSL(IT)-ABS (TSD(  IT) )  ) •GE.ZTEMP)     GO    TO    283 

KFLAG(  I T  )=  I 

TDDT( IT )=-SIGN( ADDL( IT) ,FOLUPS( IT) ) 

GO    TO    294 
263    I F(KFLAG(  IT).EQ.O  )    GO    TO    294 

IF(FOLUPS( IT)*TDDT( IT))     2  36,2  86,28  5 
Zf.5    TDDT(  IT)  =  0. 

GO    TO    2°4 
286    KFLAG( I T)=0 
294    CONTINUE 

IF(ABS( TDDT( IT) ) . GT.ADDLt  IT) )     TDDTI I T ) =S I GN ( ADDL( IT) 
*,TDDT<  IT)  ) 


IF(IMODE) 295, 297, 297 
295    CCNTINUE 

TSIL( IT)=0. 

TDCSP( IT)=0. 

TCCSK  IT)=0. 

TDCSI K IT)=TDDT( IT)/0MGTE4( IT) 
297    CCNTINUE 

TDDS(  IT)=TDDT(  IT)-FOMGTE( IT)*TSIL(  IT)-S0MGT2(  IT)* 
*    TDCSP( IT)-F0MGT3( IT)*TDCSI < IT)-0MGTE4< IT)*TOCSII( IT) 

TSIL( IT)=TSIL( IT)+TD0S( IT)*DT 

TDCSP( IT)=TDCSP(IT)+TSIL( IT)*DT 

TOCSK  IT)  =  TDCSK  IT)+TDCSP(  IT)*OT 

TDCSI I ( IT )=TDCSI I( ITJ+TDCSI (IT)*DT 
C 

C  TRANSLATIONAL  POSITION  COMMAND  OUTPUT 

C  AND  LEAD  COMPENSATION 


TSDUT  )=TDDS(  IT)*OMGQSQ(  IT)+TSIL(  IT)*GLEAD(  IT)  + 
*  TDCSP(  IT)+TIC(  IT) 

3C0  CCNTINUE 

60   CONTINUE 

RETURN 

END 
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£       llr?RSyTIye    T0     INITIALIZE    THE    FILTER    PARAMETERS    FOR 
C       THE    DRIVE    PROGRAM    (SUBROUTINE    SOIMOT) 

SLBROUTINE    SOIMIC(DTP) 

C  PROGRAMMED    BY    D.    ASTILL  (CSC)  12    SEPT     1975 

CCMMON/WSHOUT/AKRC3)  ,AK1(3)  , OMGR ( 3 ) ,OMGT  (3)  tZETR(3) 
*t     ZETT(3) 
C 

COMMON/ MOT  I ON/OMGRL( 3)  ,OMGRQ(3)  ,GMGRE(3)  ,OMGTE(3) 
It     ZETAQ(3)»ZETRL(3)rZDLF(3) ,GNRS (3),AKRL(3)  ,SOMGR( 3) 

2,  ZCMGR(3)tOMGQSQ(3) ,TOMGRE(3)  ,OMGRE2<3)  ,CMGT2(3) 

3,  TOMZT(3),FOMGTE(3) , SOMGT2<3) ,FOMGT3(3) ,0MGTE4(3) 
4»    0MGTQ(3) ,GLEAO(3) 


C 
C 


DO    200    IN-1,3 

AKRL( IN)  =  OMGRL(IN ) / ( 2 . *ZETRL(  IN  )  ) 

SCMGR( IN)=OMGR( IN)**2 

ZCMGRd  N)=2.*0MGR  (  IN)*ZETR(  IN) 

OMGQSQ( IN)=i./OMGTQ( IN)**2 

GLEAD( IN)=2.*ZETAQ( I N) /OMGTQ ( IN ) 

TCMGRE( IN)=2.*0MGRE< IN) 

0MGRE2( IN)=OMGRE( IN)**2 

0^GT2( IN)=OMGT( IN)**2 

TCMZT(IN)  =  2.*0MGT(IN)*ZETT(  IN) 

FCMGTEl IN)=4.*0MGTE( IN) 

SCMGT2( IN)=6.*0MGTE( IN)**2 

FCMGT3( IN)=4.*0MGTE( IN)**3 

CMGTEM IN)=GMGTE( IN)**4 

ZDLF( IN)=1.-EXP(-DTP*2.*ZETRL( IN)*GMGRL ( IN) ) 


C 

2C0    CONTINUE 
C 

RETURN 

END 


C   BLOCK  DATA  SUBPROGRAM  TO  INITIALIZE  THE  NON-WASHOUT 
C   FILTER  PARAMETERS 


C 


BLOCK  DATA 

CCMMON/MOTION/OMGRLi  3)  tOMGPQ( 3)  ,CMGRE(3)  ,0MGTE(3) 
I  ,  ZETAQ(3)  t  ZETP.L(3  )TZDLF(  3  )  ,  GNRS  (3),AK*L ( 3)  ,SOMGR(  3) 

2,  ZCMGR (3  ),CMGQSQ (3)  , T0MGRE13)  , 0MGRE2 ( 3 ) , CMGT2 ( 3 ) 

3,  TCMZT( 3) ,FOMGTE (3) , S0MGT2( 3) ,FCMGT3(3) ,0MGTE4<3 ) 

4,  0MGTGC3) tGLEAD(3) 

DATA  Z ETRL/ 3*. 70 7 /, ZETAQ/. 7 T. 7 ,50./ 

DATA  GNRS/ 3*1. /,0MGRL/3*1 . 2/ , OMGRQ/8.  5  ,  7. , 8.5/ 

DATA  OMGTQ/ A. 1, 7. ,  130. /,0MGTE/3*. 05/ 

CATA  0MGRE/3*.01/ 

END 
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C   CRMSBY  HUMAN  DYNAMIC  ORIENTATION  MODEL  SUBROUTINE 
C 

SUBROUTINE  HDO( XDT, NSP, SF, W , IOUT ,N ) 
C 

DIMENSION  TC<4,4)  ,TO(3,3) ,TPC( 3,3)  ,TPO<2  ,2)  ,CC(4) 
ltC0(3),XCH(4),YCH(4),ZCH(4)  ,CS(  3),XOH(3)  ,YOH(3)  ,  ZOH(3) 
2,XC(3),YC<  3),ZC(3) 7XO(2) ,GKO(3) ,GKS(3) ,ZO(3 ),CTC( 3,3) 
3,CTO<3,  3) , A(3) ,GKC(4) ,TWH(3i  , TAH(3) ,TWS<  3)  , TAS13) 
4,A0(3),  WO<  3  ),E.<JS(3  )  ,EWH(3)  ,  EAH(  3  J  ,  EAS  (  3  )  t  WNCO  (  3  ) 
5, WNCL(3)  ,VO(3) ,DVC(3) ,DVO(2),WOF(3) ,YO( 3) ,DCLD(3) 
6,WSFO(3),FN(3),Y(3) 

DIMENSION    W(N,3) f SF(N,3) 


C 

c 


CATA    DT/1.0/,NITP/10/ 

REMIND     IOUT 
REWIND    4 


NOT=NSP*OT/NITP 

lNDT=<NDT-2)/2 

WRITE ( IOUT,400)DT,INDT 
C 
C   CANAL  SPECIFICATIONS. 

DC  10  1=1,4 

READ (4, 53  0)       TC( 1,1  ),TC( I,2),TC(  I,3),TC<  I»4) 
10  CONTINUE 

DO  15  1=1,3 

READ(4,530)   TPC(  I  ,  1 ) , TPC ( I , 2 )  ,T PC ( I , 3 )  , DVC <  I  ) 
15     CONTINUE 

READ (4, 530)   CC( 1 ),CC(2) ,CC(3) ,CC<4) 

REA0(4, 530)   GKC( 1 ) ,GKC(2) ,GKC( 3) , GKC ( 4 ) 

READ<4, 550)       F SC C , TSCC , SSCC 

CALL    EULER(FSCC, TSCC, SSCC, CTC ) 
C 
C       OTOLITH    SPECIFICATIONS. 

DC  20  1=1,3 

RE AD (4, 55  0)   TO ( I  ,  1  ) , TO ( I , 2 ) , TO (  I  ,  3 ) 
20     CONTINUE 

DC  25  1=1,2 

READ(4, 550)   TPOi 1,1) ,TPO( 1,2) ,DVC(I) 
25     CONTINUE 

RE AD (4, 5  50)CO( 1 ) ,C0(2) ,C0(3) 

READ (4,  5  70)   GKO(  1) ,GK0(2) ,GKO(  3) ,GKS( 1 )  ,GKS<2)  ,GKS(3) 

READ(4,570)   FOTO , TOTO , SOTO , S ACF AC, OSPG, DFAC 

CALL  EULER(cOTO, TOTO, SOTO, CTO) 

DO  27  1=1,3 

CS(I) =CC( I )*SACFAC 
27     CONTINUE 
C 
C   INITIALIZATION 

DO  32  1=1,4 

REAC(4, 550)XCH( I ) ,YCH( I ) ,ZCH( I ) 
32     CONTINUE 

A(  1)  =  XCH(4) 

A( 2)=YCH(4) 

A(3)=ZCH(4) 

CALL  COTRN( A, CTC,  1,W0) 

DC  36  1=1,3 

READ(4,550)   XOH(  I  )  , YOH ( I )  ,  ZOH (  I ) 
36     CONTINUE 

A(  1)=X0H( 3) 

A(  2)  =  Y0H(  3) 

A(3)=ZQH<3 ) 

CALL  COTRN( A,CTG, 1 ,A0) 

DO  45  1=1,3 

RE£0(4,540)XC( I) , YC ( I ) ,ZC(I ) ,WNCO( I J,WNCL( I ) 
45     CONTINUE 

DC  50  1=1,2 

READ(4,550)   XO ( I  )  , YO ( I  )  , ZO (I ) 
50     CONTINUE 

REAC(4,550)D0LD(1  ),D0LD(2 ),D0lD(3) 

RE AD (4, 5  50) wSFOi 1 ) , WSFC ( 2 ) , WSFO ( 3 ) 
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READ(4,550)W0F(1)  ,W0F(2),W0F(3) 

READ<4,  5  80)    TDVEL,T0POS,TNC,FNOISE 

FP=1.0-EXP( -DT/T0POS) 

FC=TDVEL 

FN<  1)  =  EXP(-DT/T,NC) 

FN<2)=T,\K><  l.-FN(  1  )  )/OT-FN(  1) 
55  FN(3)  =  1.-TNC*(1.-*N<  1))/DT 

C 

C       MAIN    PROGRAM    CYCLE. 
C 

C  CCMPUTE    CURRENT    STIMULUS    IN    HEAO    COORDINATES 

C  (EVERY    DT/NITP    SEC.) 

C  1.    ANGULAR    ROTATION    VECTOR     (TWHJ     AT    (TIME). 

C  2.    SPECIFIC    FORCE    VECTOR    ( TAH )     AT    (TIME+OT/2) 

C  3.    TRUE    DOWN    VECTOR-   AT    (TIME+DT/2). 


C 

C 
C 


DO    450     ITIME-1, I  NOT 

DO    100    I=1,NITP 

TIME=( ITIME-1 >*DT+ I *DT/NITP 

J=2C.*TIME+1.5 
K=J+10 


TAH(  1)  =  SF( K,i)/32.2 

TAH(2)=SF(K,2)/32  .2 

TAH(3)=SF(K,3)/32.2 

TWHU  )  =  W(  J,  1) 

TWH(2)=W( J, 2) 

TWH(3)=W( J, 3) 
C 
C     TRANSFORM  TO  SENSOR  COORDINATES. 

CALL  COTRN(TWH,CTC,0,TWS) 

CALL  COTRN(TAH,CTO,0, TAS) 
C 

C     SENSOR  STIMULATION  (EVERY  DT/NITP  SEC): 
C       LSING  CURRENT  STIMULUS  VALUES,  UPDATE  STATE  VECTORS 
C       FOR  3  CANALS  IXC,  YC,  AND  ZC),  AND  3  OTOLITHS 
C       <X0,  YOt  ZO),  ANO  COMPUTE  AFFERENT  FIRING  RATES 
C       (CSX,  CSY,  CSZ). 
C 

S=TWS(1 ) 

CALL  SVUP0( XC,TPC,DVC,S,CSX,CC,3,4) 

S  =  TWS( 2  ) 

CALL  SVUPD( YC,TPC,DVC,S,CSY,CC, 3,4) 

S=TWS(3) 

CALL  SVUPD( ZC,TPC,DVC,S,CSZ,CC,3,4) 

S=TAS( 1 ) 

CALL  SVUPD( XO,TPO,DVOtS,CSX,CO, 2,3) 

S=TAS(2) 

CALL  SVUPD( YO,TPO,DVC,S,OSY,CO,2,3) 

S  =  TAS(  3  ) 
100  CALL  SVUPD( ZO ,TPO , DVG, S ,OSZ , CS , 2,3) 
C 

C     OPTIMAL  ESTIMATOR  (UPDATE  EVERY  DT  SEC): 
C       GET  CANAL  AND  OTOLITH  SYSTEM,  STATE  ESTIMATES  FROM 
C       STEADY  STATE  KALMAN  FILTERS. 
C 

CALL  SSKF(XCH,CSX,TC,CC,GKC,4) 

CALL  SSKF( YCH,CSY,TC,CC,GKC ,4) 

CALL  SSKF(  ZCH,CSZ  ,TC,CCGKC,4) 

CALL  SSKP(X0H,0SX,T0,CO,GKO,3) 

CALL  SSKF( Y0H,GSY,T0,CC,GK0,3) 

CALL    SSKF(  ZGH,OSZ  ,T0, CS , GKS , 3 ) 
C 
C  ENTER    ROTATION    RATE    ESTIMATE    VECTOR    (CANAL    ESTIMATE). 

EVS( 1)=XCH(4) 

EkS(2)  =  YCH<4) 

E'*S(3  )^ZCH(4) 
C 
C  ENTER    SPECIFIC    FORCE    ESTIMATE    VECTOR    (OTOLITH    ESTIMATE) 

EAS(1)=X0H(3) 

EAS( 2)=YOH(3) 


84 


c 

C     SACCULE  NON-LINEARITY.   - 

EAS(3)=AMAX1<  .6*  (  ZOH(  3  )  +.41  69  )-  .4169,  -  .  4  169  ) 

9  §!lT^£E^A5fyIIUE0    0F    OTOLITH    ESTIMATE    TO    VALUE    HELD 

£  BEFORE    CONSIDERATION    OF    SACCULE    NON-LINEARITY. 

C  (THEREFORE,     NON-LINEARITY    EFFECTS    ONLY    DIRECTION    OF 

C  THE    OTOLITH    ESTIMATE). 

C 

CALL  NORM(EAS,Y) 

DC  130  1=1,3 
130  EAS( I)=Y( I)*SQRT(XOH( 3)**2+YOH< 3)**2+Z0H<3)**2) 

C     TRANSFORM  TO  HEAD  COORDINATES 

CALL  COTRN( EWS,CTC,l ,EWH) 

CALL  COTRN( EAS , CTO, 1 , E AH) 
C 

C  DOWN  AND  W  ESTIMATOR  (UPDATE  EVERY  DT  SECONDS). 
C  CCMBINE  OTOLITH  AND  CANAL  ESTIMATES  TO  FORM  NEW 
C     ESTIMATE  OF: 

C       1.  PERCEIVED  DOWN  (DNEW)  AT  (TIME+DT/2). 
C       2.  PERCEIVED  ACCELERATION  (ACC)  AT  (TIME+DT/2). 
C       3.  PERCEIVED  ANGULAR  VELOCITY  ( WTOT )  AT  (TIME). 
C 

44C  CALL  DCWN(DOLD,E*'H,EAH,AO,WSFO,FD,  CT  ,TDPOS  ,  DFAC  ,  WOF 
*,  WNCO ,WNCL ,FN,IOUT) 

45C  CCNTINUE 
C 
C   FORMAT  STATEMENTS. 


C 


400  FCRMAT( F6.3, 14) 
52C  FCRMAT*  E15.8, 15 ) 
5  50  FCRMAT( 4E15.3) 
540  FCRMAT( 5E15.8) 
550  FORMAT( 3E15.8) 
57C  FORMAT( 6E 1 2.5) 
580  FCPMAT( 4E15.3  ) 

RETURN 
END 


C   DCWN  ESTIMATOR  AND  W  ESTIMATOR 
C 

SUBROUTINE  DOWN( DOLD,WN,SN,SO,WSFO,T,DT, TDPS , DFAC , WOF 
*f  WNCO  ,WNL,FN, I  OUT) 
C 

C   COWN  IS  DETERMINED  SY  RELYING  CN  LOW  FREQUENCY  OTOLITH 
C   ESTIMATES,  CANAL  ESTIMATES  THAT  ARE  CONSISTENT  WITH 
C   HIGH  FREQUENCY  OTOLITH  ESTIMATES,  AND  THE  HIGH  FREQUENCY 
C   PCRTICN  OF  CANAL  ESTIMATES  NOT  CONFIRMED  6Y  OTOLITHS. 
C 
C 
C 

c 
c 
c 
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F(2)=TD 
F(3)=l. 

CALL  CR 
CALL  NO 
CALL  VA 
DC  10  I 
WSF(I)= 
WCF( I) = 
WCD(  I  )  = 

10  WSFO(  I  ) 
WODM=SQ 
CALL  NO 
WCPWO=W 
IF( WCPW 

11  WCPWD=0 

12  WMAG=-W 
DC  13  I 
X(  I)=WO 
WNC< I )= 
WNL( I)= 
WNCC(  I) 
WNCH(  I) 

13  X(  I  )  =  X( 
CALL  CR 
CALL  NC 
WCPARM= 
DC  15  I 
ROTO( I ) 

15  HRCTOU 
CALL  RO 
WCPARM= 
CC  20  I 
RSCC(  I) 

20  RTOT(  I  ) 
CALL  RC 
C0LL  VA 
PHI=FPO 
CALL  CR 
CALL  NO 
DC    30     I 

30  RPGSi I) 
CALL  RC 
CALL  NO 
DO    40     I 

40  X(  I)=DO 
CALL  NO 
WPARM=W 
CALL  CR 
CALL  NO 
CALL  VA 
DC  50  I 
WPERP(  I 
WPAREU 
DCLC( I) 
SC(  I)=S 
UTOT( I ) 
ACC(  I)  = 

50    CCNTINU 

45    WRITE<  I 

*, WT0T<3 

ICC    FCPVAK 

RETURN 
END 


VEL*(1.-F( 

-T0VEL*(1. 

OSS(SO,SN, 

RM( WSF,X) 

NG( SO,SN,A 

=  1,3 

ANGSF*X(  I  ) 

F(  1 )*WOF(  I 

WSF(  n-WOF 

=WSF(I) 

RT( WGD( 1 )* 

RM( WOO, WOO 

N( 1)*W0DN( 

D)     12,12,1 

.0 

CPWO*OT 

=  1,3 

DN( I  )*AMIN 

-WN(  I  J*OT- 

FN(  1  )*WNL< 

=WNC(  I  ) 

=  WNC( I)-WN 

I  )+WNCH( I  ) 

CSS     (DOLDT 

RM( ^OTO,FJ 

WOF(  1)*F<1 

=  1,3 

=WOPARM*Fl 

)=ROTO( I)/ 

TATE(DOLD, 

X(  1)*F< 1>  + 

=  1,3 

=  (  X(  D-WCP 

=RSCC( I)+R 

TATE(OOLD, 

NG( DNEW,SN 

S*FE  E 

OSS(DNEW,S 

RM(RPOS,X) 

=  1,3 

=PHI*X( I ) 

TATE(DN£W, 

RM(X,DNEW) 

=  1,3 

LD( I  1+DNEW 

RM(X,DAVG) 

N( 1)*0AVG( 

OSS(DNEW,D 

RM( WPERP ,X 

NG(DCLD,DN 

=  1,3 

)  =  x(  n-PHi 

)=WPARM*DA 

=ONEW(I) 

N(  I  ) 

=WPERP<I )+ 

DFAC*DNEW( 

E 

OUT,  100) AC 

),DNEW( 1  )  , 

9E14.7) 


1) )/DT-F(l) 
-F(l  )  )/DT 

WSF) 

NGSF) 


)+F(2)*WSF0( I)+F(3)*WSF< I) 
(I) 

WOD( 1)+W0D( 2)*WOD(2)+WOO<3 )*WOH< 3) ) 

N) 

1 )+WN(2)*W00N( 2) +WN(3)*WCDN(3) 

1 


KWMAGtWODM) 
X(  I) 

I  )  +  FN< 2)*WNC0<  I)+FN( 3)*WNC (I  ) 

HI) 

SN,ROTO) 

)+WOF(2)*F(2)+WOF( 3)*F(  3) 

I  ) 
2. 

HROTO,F) 
X(2)*F(2)+X(3)*F<  3) 

ARN!*F<  I)  ) 
0T0( I ) 
RTOT,DNEW) 
,FEE) 

N,RPOS) 


RPOS,X) 


(I) 

1 )+WN( 2)*DAVG< 2)+WN(3)*DAVG( 3) 
OLD, WPERP) 
) 
EW,PHI  ) 

/DT 
VG(  I  ) 


WPARE( I) 
II-SNU) 

C(l)  ,ACC(2 ) ,ACC(3),WT0T( 1  )  ,WTOT{ 2) 
DNEW( 25 ,QNEW( 3) 
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£       IVigEyiJNE    TO    DEVELOP    A    COST    ALGORITHM    BASED    ON    THE 
£       9l£E^E£CE    BETWEEN    THE    PERCEPTION    OF    SIMULATOR    AND 
C       AIRCRAFT    MOTION 


C 

c 
c 


SLBROUTINE    COST ( CAD, CAVG»CRMS ) 

CCMM0N/FLIGHT/XA(4384,3),XR(4884,3),P,C,  I  PAR 

D I  PENS  ION  AA(3) , AS(3),WA(3) ,WS( 3),DWNA( 3  ),DWNS( 3) 
DMA  WTH,ATH/. 0146,  .005/ 

CAVG=0. 

CRMS=0. 

CAC=0. 
C 

CV<AV  =  0. 

C4AV=0. 

CCAV=0. 
C 

CWRMS=0. 

CARMS=0. 

CCRMS=0. 
C 

ACW=0. 

ACA=0. 

ACC=0. 
C 

C       READ     IN    THE    DATA    FROM    THE    TWO    FILES 
C 

REWIND    9 

REWIND    8 
C 

RE40(9,30)0TS,NSS 

READ(8,30)DTA,NSA 
C 
C       CCMPLTATION    OF    THE    COST    ESTIMATION    VALUES 


C 


c 

c 

c 
c 


DO    40    1=1, NSA 

READ(9,10JAS(1 ) , AS ( 2 )  , AS ( 3 ) , WS (  1 ) , WS ( 2 ) , WS ( 3 )  ,DWNS(  1) 
*,DWNS( 2) ,DWNS( 3) 

REAC(8,  10) A  All ) , AA(2 )  ,AA( 3) , WA (  1  )  ,WA(2)  ,WA( 3)  ,DWNA(  1  ) 
*,DV<NA(2)  ,DWNA(3) 

CW1=WA(  1  )-WS(l  ) 
DV»2  =  WA(  2)-WS(2) 
DW3=WA( 3)-WS(3) 

D/U=AA(  1)-AS(  1  ) 
C/^2  =  AA(  2)-AS(2) 
C43=AA( 3)-AS(3) 

C01=DWNA<1 )-DWNS(  1  ) 
CC2=DWNA(2)-DWNS( 2) 
CD2=0WNA(3)-DWNS( 3) 

Dk=SQRT (DW1*DW1+QW2*DW2+DW3*DW3 ) 
CA=SQRT (DA1*DA1+DA2*DA2+DA3*DA3 ) 
0D=SQRT(DD1*DD1+DD2*DD2+DD3*DD3) 

AWA=SQRT(WA(1)*WA(1)+WA(2 )*WA( 2 ) +WA ( 3 ) *W A (  3  )  ) 
AAA=SQRT(AA(1)*AA(1)+AA(2)*AA(2)+AA(3)*AA(3)) 

Cw=DW/AMAXl (AWA, MTH) 
CA  =  CA/AMAXKAAA,  ATH) 
CC  =  CD 

CWAV=CWAV+CW 
CAAV=CAAV+CA 
CCAV=CDAV+CD 

CVsRMS=CWRMS+Crl*CW 
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c 

40 

c 


c 

70 
50 

60 

30 
10 

c 


CARMS=CARMS+CA*CA 
CCRMS=CDRMS+CD*CD 

ACW=ADW+DW 
ACA=ADA+DA 
ACC=ADD+D0 

CONTINUE 

OAV=CWAV/NSA 
CAAV=CAAV/NSA 
CCAV=COAV/NSA 
C£VG=CWAV+CAAV+CDAV 

CWRMS«SQRT(CWRMS/NSA) 

CARMS=SGRT(CARMS/NSA) 
CDRMS=SQRT(CDRMS/NSA) 
CRWS=CWRMS+CARMS+CDRMS 

ACw=AOW/NSA 
ACA=ACA/NSA 
ADO=AOD/NSA 
C*D=ADW+ADA+A0D 

WRITE<6,70)     ADW,ADAfAOO,CAD 

FCRMAT(//»  AO^SFll.6,1 

:,F11.6,<  TOTAL=' ,F10.6) 

fcRITE(6,50)CWAV,CAAV,CDAV,CAVG 

FCRMAT{/«  CWAV=»,F10.6t '       CAAV=»,F10 

■■,     F10.6,'  TOTAL  =  «  ,F10.6) 

WRITE (6 ,60) CWRMSt CARMStCDRMS,CRMS 

FCRMAT(/«  CWR^S='  ,F9.6, •        CARMS='»F9 

:,F9.6, •  TCTAL=«  ,c10.6/) 

FCR,MAT(  F6.3,I4) 

FCPMAT(9E14.7) 

RETURN 

END 


A0A*'vF11.6f<       ADD  =  « 


6t 


6, 


CDAV=* 


CDRMS=I 
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C   SUBRCITINE  TO  DIFFERENTIATE  THE  TIME  SERIES 
C   CONSISTING  OF  THE  ROTATIONAL  POSITION 
C   COMMANDS  FROM  THE  DRIVE  PROGRAM 


C 


SUBROUTINE    DI FF ( DT , Y, YD, NSP,N ) 


(DIMENSION    Y(N,3)  ,YD(N,3) 

DC    10    1=1, NSP 

DO    40    J  =  l,3 

I F  t 1-1 ) 20, 20, 50 
50  IF(I-NSP) 60,30,30 

20  YC<  I,J)=<Y<  1  +  1, J )-Y(  I,J))/DT 

GC    TO    40 
30  YCU,  J)=(  Y(  I,J)-Y(I-1,J))/DT 

GC    TO    40 
60  VC( I, J)=(Y(  1  +  1 ,J)-Y(I-1,J) ) /(2.*DT) 

40  CONTINUE 

10  CONTINUE 

C 

RETURN 

ENC 


C       SUBROUTINE    TO    TRANSFORM    THE    COORDINATE    SYSTEM    OF    A    VECTOR 
C 

SUBROUTINE    COTRN ( A , B, N ,C ) 
C 

C       COORDINATE    TRANSFORM: 
C 

C  FRCM    HEAD    TO    SENSOR    IF    N=0     (DIRECT    TRANSFORMATION) 

C  FROM    SENSOR     TO    HEAD    IF    N=-l     (INVERSE    TRANSFORMATION) 

C  FRCM    SENSOR    TO    HEAD    IF    N=l     (INVERSE    TRANSFORMATION 

C  LSING    THE    TRANSPOSE    FOR    ORTHOGONAL    MATRICES) 

C 

C  A=CRIGINAL    VECTOR 

C  B=TRANSFORMATION    MATRIX 

C  C=TRANSFORMED    VECTOR 

C 

DIME NS I  ON  A(3),B(3,3),C(3),BI(3,3),L(3),M(3) 
C 

IF(N)20,10,30 
C 
10     DC  15  1=1,3 

C( I )=B(  I, 1) *A(  1 )  +  6< I , 2)*A(2 )  +  B( I  ,  3)*A(3) 
15     CONTINUE 

RETURN 
C 
20     CALL  MCPY(3,6I,3, 3,0) 

CALL  MINV(BI  ,3 ,D, L,M) 

DO  35  1=1,3 

C(I )=BI  ( 1,1 >*A( 1  ) +311  I,2)*A( 2)+BI (I, 3)*A( ?) 
35     CONTINUE 

RETURN 
C 
30     DC  40  1=1,3 

C (I)  =  B(  1, I  )*A(1 )+B(2,  I J*A(2  )+B(  3,1  )*A( 3  ) 
40     CCNTINUE 
C 

RETURN 

END 
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SUBROUTINE  EULER ( F ,T, S , CT ) 

C   PRODUCE  DIRECTION  COSINE  MATRIX  (CT)  GIVEN  EULER  ANGLES 

C   (F,T,  AND  S) 

C 

DIMENSION  CT(3t3) 

CT(  l,l)=COS(S)*COS(F)-COS(T)*SIN(F)*SIN(S) 

CT(2,1)=-(S IN(S)*COS(F)+COS(T)*SIN(F)*COS(S)) 

CT(3,1 )=SIN(T)*SIN(F) 

CT(l,2)=COS(S)*SIN(F  )+COS ( T } *COS < F )*S IN ( S ) 

CT(2,2)=C0S(T )*COS(F)*COS(S)-SIN(SJ*SIN<F) 

CT(3,2)=-SIN(T)*COS(F) 

CT( 1,3)=SIN(T)*SIN(S) 

CT(2,3)=SIN(T)*COS(S) 

CT(3,3)=C0S(T) 
C 

RETURN 

END 


SLBROUTINE  CROSS  (A,B,C) 

C 

C    C  =  A  X  B 

C 

DIMENSION  A(3> ,B( 3) ,C( 3) 
C(1)=A( 2)*B(3)-A( 3)*B( 2) 
C(2)  =  A( 3)*B(1)-A<  1)*9( 3) 
C(3)=A<  1)*B(2)-A(  2)*B(  1) 

C 

RETURN 
ENC 


SUBROUTINE  VANGC Af B, PHI ) 

C 

C       °HI     =    ANGLE    BETWEEN     A    AND    B. 


C 


DIMENSION    A(3 ) t 3(3), AN( 3) ,3N( 3) 

CALL    NOR*1(  A, AN) 

CALL    NORM( 3,BN) 

PH£*ARCOS (AN(1 )*BN(1 )+AN( 2 ) *BN ( 2 ) +AN ( 3 ) *BN<  3) ) 

RETURN 
ENC 


SUBROUTINE  NORM(A,AN) 

C   AN  =  UNIT  VECTOR  IN  DIRECTION  OF  VECTOR  A. 

C 

DIMENSION  A(3),AN(3)     ,      ,.,   ,_%1 
AV  =  SCRT(A(1)*A(1)+A(2)*A(2)+A(;>)*A(:>)) 

IFUM-l.E-06)     5t5,6 

5  Af=l.E-06 

6  DC     10    1=1 i 3 
1C    AN(I)=A  (  D/AM 

C 

RETURN 
END 


90 


SUBROUTINE  SVUPD ( XtT , D, S, Y ,C , Nt M ) 

C  STATE  VECTOR  UPDATE: 

C 

C  X(NEW)  =T*X(OLD)  +  D*S 

C  Y(NEW)  =  C*X(NEW) 

C  WHERE:  N  IS  DIMENSION  OF  STATE  VECTOR 

C  X  IS  STATE  VECTOR 

C  T  IS  TRANSITION  MATRIX 

C  D  IS  DRIVING  VECTOR 

C  S  IS  STIMULUS 

C  Y  IS  OUTPUT  (AFFERENT  FIRING  RATE) 


C 


DIMENSION    X(N),T(N,N),D(N),C(M),R(9) 
DC    5    1=  1  • N 
5    R(  I)  =  X(  I) 
DC     10    1=1, N 

xm-o<  n*s 

DC  10  J=1,N 
10    X(  I)=X(  I)+T{ I,J)*R(J) 

Y=C(M)*S 

CC  20  I = 1 , N 
20    Y=\+C( I )*X(  I) 

RETURN 
END 


SLBROUTINE    SSKF  (  XH  ,  Y  ,  TM  tC  ,  GK,  N  ) 
C 

C  STEACY    STATE    KAlMAN    FILTER    (UPDATE    EVERY    DT    SECONDS) 
C 

C  WHERE:     XH    IS    STATE    VECTOR    ESTIMATE 

C  TM    IS    TRANSITION    MATRIX 

C  GK     IS    KALMAN    GAIN    MATRIX 

C  Y    IS    SENSOR    SYSTEM    OUTPUT 

C  CIS    OUTPUT    MATRIX 


C 


DIMENSION    XH(N),TM(N,N),C(N),GK(N),S(9) 

DO    40    1=1, N 

S(  I) =0.0 

DC  40  J=1,N 
40    S (  I  )  =  S (  I)+TM( I, J)*XH( J) 

EV=0.0 

DC  45  1=1, N 
45    EN=EM+S ( I )*C( I ) 

DC  50  J=1,N 
50    XH(  J)  =  S(  J)+GK(J)*(Y-EM) 

RETURN 
ENC 


SLBROUTINE    ROTATE     (A,R,AR) 
C 

C       AR    =     A    ROTATED    ABOUT    R    BY    ANGLE    (RAD)    EQUAL    TO    THE 
C       MAGNITUDE    OF    R 


C 


CINENS ION  A(3),R<3),AR(3),AP(3),APN<?) 
CALL  CROSS  (R,A,AP) 
CALL  NORM  (APTAPN) 

AMAG=SQRT( A(l )*a(  1  )  +A  (  2  5  *A  (  2  )  +A  (  3  )  *A  (  3  )  ) 
PHI=SQRT(R( 1)*R( 1 )+R(2)*R(2)+R( 3)*R(3) ) 
DC  10  1=1,3 
10  AR ( I)=AMAG*SIN( PHI )*APN(I ) +COS (PHI)*A{ I ) 

RETURN 

ENC 
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